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ABSTRACT 

Various different physical processes contribute to the star formation and stellar mass assembly histories of galaxies. One important 
approach to understand the significance of these different processes on galaxy evolution is the study of the stellar population content of 
today’s galaxies in a spatially resolved manner. The aim of this paper is to characterize in detail the radial structure of stellar population 
properties of galaxies in the nearby universe, based on a uniquely large galaxy sample considering the quality and coverage of the data. 

The sample under study was drawn from the CALIFA survey and contains 300 galaxies observed with integral field spectroscopy. 

These cover a wide range of Hubble types, front spheroids to spiral galaxies, while stellar masses range from M* ~ 10 9 to 7 x 10 11 
M q . We apply the fossil record method based on spectral synthesis techniques to recover the following physical properties for each 
spatial resolution element in our target galaxies: the stellar mass surface density (//*), stellar extinction ( A v ), light-weighted and 
mass-weighted ages ((log age) L , (log age) M )* and mass-weighted metallicity ((logZ*) M ). To study mean trends with overall galaxy 
properties, the individual radial profiles are stacked in seven bins of galaxy morphology (E, SO, Sa, Sb, Sbc, Sc and Sd). We confirm 
that more massive galaxies are more compact, older, more metal rich, and less reddened by dust. Additionally, we find that these trends 
are preserved spatially with the radial distance to the nucleus. Deviations from these relations appear correlated with Hubble type: 
earlier types are more compact, older, and more metal rich for a given M*, which evidences that quenching is related to morphology, 
but not driven by mass. Negative gradients of (log age) L are consistent with an inside-out growth of galaxies, with the largest (log age)i 
gradients in Sb-Sbc galaxies. Further, the mean stellar ages of disks and bulges are correlated, with disks covering a wider range of 
ages, and late type spirals hosting younger disks. However, age gradients are only mildly negative or flat beyond R ~ 2 HLR, indicating 
that star formation is more uniformly distributed or that stellar migration is important at these distances. The gradients in stellar mass 
surface density depend mostly on stellar mass, in the sense that more massive galaxies are more centrally concentrated. Whatever 
sets the concentration indices of galaxies obviously depends less on quenching / morphology than on the depth of the potential well. 

There is a secondary correlation in the sense that at the same M* early type galaxies have steeper gradients. The /r* gradients outside 
1 HLR show no dependence on Hubble type. We find mildly negative (logZ*) M gradients, shallower than predicted from models of 
galaxy evolution in isolation. In general, metallicity gradients depend on stellar mass, and less on morphology, hinting that metallicity 
is affected by the depth of both - potential well and morphology/quenching. Thus, the largest (logZ*) M gradients occur in Milky 
Way-like Sb-Sbc galaxies, and are similar to those measured above the Galactic disk. Sc spirals show flatter (logZ*) M gradients, 
possibly indicating a larger contribution front secular evolution in disks. The galaxies from the sample have decreasing-outwards 
stellar extinction; all spirals show similar radial profiles, independent from the stellar mass, but redder than E's and SO’s. Overall we 
conclude that quenching processes act in manners that are independent of mass, while metallicity and galaxy structure are influenced 
by mass-dependent processes. 

Key words. Techniques: Integral Field Spectroscopy; galaxies: evolution; galaxies: stellar content; galaxies: structure; galaxies: 
fundamental parameters; galaxies: bulges; galaxies: spiral 


1. Introduction 

Galaxies are a complex mix of stars, gas, dust, and dark matter, 
distributed in different components (bulge, disk, and halo) whose 
present day structure and dynamics are intimately linked to their 
assembly and evolution over the history of the Universe. Differ¬ 
ent observational and theoretical approaches can be followed to 
learn how galaxies form and evolve. 

Theoretically, the formation of large-scale structures arise 
through the evolution of cold dark matter. In this picture, small- 


scale density perturbations in the dark matter collapse and form 
the first generation of dark matter halos, that subsequently merge 
to form larger structures such as clusters and superclusters 
(Springel et al. 2005; De Lucia et al. 2006). This basic hierarchi¬ 
cal picture is able to explain the global evolution of the star for¬ 
mation rate density of the universe, with galaxy peak formation 
epoch at redshift 2-3 (e.g. Madau & Dickinson 2014, and refer¬ 
ences therein). The stellar components formed at earlier epochs 
likely evolved into elliptical galaxies and bulges through merg¬ 
ers of the primordial star-forming disks (Elmegreen et al. 2007; 


Article number, page 1 of 42 




A&A proofs: manuscript no. Califa5_metallicity_vl2_8Jun 


Bournaud et al. 2007). However, this framework fails to explain 
how the galaxy population emerges at z ~ 1, and how the present 
day Hubble sequence of galaxies was assembled. 

The growth of galaxies is not related in a simple way to the 
build up of dark matter; the interplay of energy and matter ex¬ 
change (between the process of gas accretion and cooling and 
star formation) is essential to grow the gaseous and stellar com¬ 
ponents in galaxies. Feedback mechanisms resulting from stel¬ 
lar winds, supernova explosions, and AGN are relevant to stop 
the gas collapse and cooling, quenching the star formation and 
hence galaxy growth (Silk & Rees 1998; Hopkins et al. 2011). 
Although these processes are difficult to implement in theoreti¬ 
cal models, they are essential to explain the masses, structures, 
morphologies, stellar populations, and chemical compositions of 
galaxies, and the evolution of these properties with cosmic time. 

Recently, a new set of cosmological hydrodynamic simula¬ 
tions have started to predict how the spatially resolved informa¬ 
tion of the properties of stellar populations in galaxies can con¬ 
strain the complex interplay between gas infall, outflows, stellar 
migration, radial gas flows, and star formation efficiency, in driv¬ 
ing the inside-out growth of galactic disks (Brook et al. 2012; 
Gibson et al. 2013; Few et al. 2012; Pilkington et al. 2012a; 
Minchev et al. 2014). Radiative cooling, star formation, feed¬ 
back from supernovae, and chemical enrichment are also in¬ 
cluded in simulations to predict radial metallicity gradients as 
a function of merging history. Shallow metallicity gradients are 
expected if elliptical galaxies result from major mergers (e.g. 
Kobayashi 2004), but a minor merger picture for the formation 
of ellipticals can successfully explain the strong size evolution of 
massive galaxies (Oser et al. 2012). This late-time accretion of 
low mass and metal poor galaxies (dry mergers) into the already 
formed massive galaxy can produce a variation of the age and 
metallicity radial structure of the galaxy as it increases in size. 
Galactic stellar winds and metal cooling have also an important 
effect on these ex-situ star formation models, predicting different 
behaviour of the mass and metallicity assembly in massive early 
type galaxies, and in the radial gradient of present stellar popu¬ 
lations properties of galaxies (Hirschmann et al. 2013, 2015). 

In summary, these theoretical works show that observational 
data with spatial information of the mass and metallicity assem¬ 
bly and their cosmic evolution, and the present radial structure 
of stellar population properties (stellar mass surface density, age, 
metallicity) contain relevant information to constrain the forma¬ 
tion history of galaxies, and the physics of feedback mechanisms 
involved. 

Observationally, a first step is to study what kinds of galax¬ 
ies are there in the Universe, and which are their physical prop¬ 
erties. Attending to their form and structure, galaxies can be 
grouped into a few categories. Results show that most of the 
massive nearby galaxies are ellipticals, SO’s, or spirals (Blanton 
& Moustakas 2009) following the Hubble tuning fork diagram. 
In this scheme, SO’s are a transition between spirals and ellipti¬ 
cals (Cappellari et al. 2013), and the bulge/disk ratio increases 
from late to early type spirals. At the same time, galaxy prop¬ 
erties such as color, mass, surface brightness, luminosity, and 
gas fraction are correlated with Hubble type (Roberts & Haynes 
1994), suggesting that the Hubble sequence somehow reflects 
possible paths for galaxy formation and evolution. However, the 
processes structuring galaxies along the Hubble sequence are 
still poorly understood. 

Integral Field Spectroscopy (IFS) enables a leap forward, 
providing 3D information (2D spatial + ID spectral) on galax¬ 
ies. Such datacubes allow one to recover two-dimensional maps 
of stellar mass surface density, stellar ages, metallicities, extinc¬ 


tion and kinematics, as well as a suit of nebular properties such 
as gas kinematics, metallicity, excitation, and etc. Until a few 
years ago IFS was used to target small samples of galaxies. De¬ 
tailed programs such as SAURON (Bacon et al. 2001), VENGA 
(Blanc et al. 2013), (U)LIRs at z <0.26 (Arribas et al. 2010), 
PINGS (Rosales-Ortega et al. 2010), or DiskMass Survey (Ber- 
shady et al. 2010), have been limited to less than a hundred 
galaxies, but it is more than fair to recognize that to get these 
amounts of IFU data was a challenge at the time. ATLAS 3D 
(Cappellari et al. 2011) represented a step forward, with the ob¬ 
servation of a volume-limited sample of 260 galaxies, but with 
three important limitations; the sample only includes early-type 
galaxies, the field of view is limited to 1 effective radius, and the 
spectral range is restricted from H/i to [NI]T5200. 

CALIFA (Calar Alto Legacy Integral Field Area) is our on¬ 
going survey of 600 nearby galaxies at the 3.5m at Calar Alto 
(Sanchez et al. 2012) 1 . The data set provided by the survey (see 
Husemann et al. 2013 for DR1; Garcfa-Benito et al. 2015 for 
DR2) is unique to advance in these issues not only because of 
its ability to provide spectral and spatial information, but also 
because: a) It includes a large homogeneous sample of galaxies 
across the color-magnitude diagram, covering a large range of 
masses (10 9 to 10 12 M 0 , Gonzalez Delgado et al. 2014c), and 
morphologies from Ellipticals (E0-E7), Lenticulars (SO-SOa), to 
Spirals (Sa to Sm) (see Walcher et al. (2014) for a general de¬ 
scription of the sample), b) It has a large field of view (74"x65") 
with a final spatial sampling of 1 arcsec, and a resolution of 
~ 2.5 arcsec, allowing to spatially resolve well the stellar pop¬ 
ulation properties, and to obtain the total integrated properties, 
such as galaxy stellar mass, and stellar metallicity. c) It covers 
the whole rest-frame optical wavelength at intermediate spectral 
resolution, including the most relevant absorption diagnostics for 
deriving the stellar population properties. 

Previous papers in this series have used the first -100 dat¬ 
acubes of the survey to derive spatially resolved stellar popula¬ 
tion properties by means of full spectral fitting techniques. We 
have obtained that: 1) Massive galaxies grow their stellar mass 
inside-out. The signal of downsizing is shown to be spatially pre¬ 
served, with both inner and outer regions growing faster for more 
massive galaxies. The relative growth rate of the spheroidal com¬ 
ponent (nucleus and inner galaxy), which peaked 5-7 Gyr ago, 
shows a maximum at a critical stellar mass M* ~ 7 X 10 10 M o 
(Perez et al. 2013). 2) The inside-out scenario is also supported 
by the negative radial gradients of the stellar population ages 
(Gonzalez Delgado et al. 2014c). 3) Global and local relations 
between stellar mass, stellar mass surface density and stellar 
metallicity relation were investigated, along with their evolu¬ 
tion (as derived from our fossil record analysis). In disks, the 
stellar mass surface density regulates the ages and the metallic¬ 
ity. In spheroids, the galaxy stellar mass dominates the physics 
of star formation and chemical enrichment (Gonzalez Delgado 
et al. 2014c,a). 4) In terms of integrated versus spatially resolved 
properties, the stellar population properties are well represented 
by their values at 1 HLR (Gonzalez Delgado et al. 2014c,a). The 
CALIFA collaboration has also compared the age and metallicity 
gradients in a subsample of 62 face-on spirals and it was found 
that there is no difference between the stellar population prop¬ 
erties in barred and unbarred galaxies (Sanchez-Blazquez et al. 
2014). 

In this paper we extend our study of the spatially resolved 
star formation history of CALIFA galaxies to derive the radial 
structure of the stellar population properties as a function of 


1 http://califa.caha.es 
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Hubble type, and galaxy stellar mass, M*. The goals are: 1) 
To characterize in detail the radial structure of stellar population 
properties of galaxies in the local universe. 2) To find out how 
these properties are correlated with Hubble type, and if the Hub¬ 
ble sequence is a scheme to organize galaxies by mass and age, 
and/or mass and metallicity. 3) To establish observational con¬ 
straints to galaxy formation models via the radial distributions 
and gradients of stellar populations for disk and bulge dominated 
galaxies. 

This paper is organized as follows: Section 2 describes the 
observations and summarizes the properties of the CALIFA 
galaxies analyzed here. In Sec. 3 we summarize our method for 
extracting the SFH, based on the fossil record method, and we 
explain the main differences between the analysis presented here 
and that in previous works. Sec. 4 presents results on the galaxy 
stellar mass, half light and half mass radii (HLR, HMR, respec¬ 
tively), and galaxy averaged stellar metallicity. Sec. 5 deals with 
the spatially resolved properties of the stellar population: stel¬ 
lar mass surface density, /r*; luminosity weighted mean age, 
(log age)/ \ mass weighted mean metallicity, (log Z*) M ; and stel¬ 
lar extinction. Ay. We discuss the results in Sec. 6; and Sec. 7 
presents the conclusions. 


2. Sample, and Observations, data reduction 

2.1. Sample and morphological classification 

The CALIFA mother sample consists of 939 galaxies selected 
from SDSS survey in the redshift range z = 0.005-0.03, and 
with r-band angular isophotal diameter of 45-80". These criteria 
guarantee that the objects fill well the 74"x64" FoV. The sample 
includes a significant number of galaxies in different bins in the 
color-magnitude diagram (CMD), ensuring that CALIFA spans 
a wide and representative range of galaxy types. 

The galaxies were morphologically classified as Ellipticals 
(EO-7), Spirals (SO, SOa, Sab, Sb, Sbc, Sc, Scd, Sd, Sm), and 
Irregulars (I). The classification was carried out through visual 
inspection of the r-band images averaging the results (after clip¬ 
ping outliers) from five members of the collaboration. Galaxies 
are also classified as B for barred, otherwise A, or AB if it is 
unsure, and as M if it shows "merger" or "interaction features" 
(Walcher et al. 2014). 

The sample for this paper comprises the 312 CALIFA galax¬ 
ies observed in both V1200 and V500 setups as of January 2014. 
The 12 galaxies showing "merger or interacting features" are not 
discussed here, leaving a main sample of 300 objects with a well 
defined morphology. For this work we have grouped galaxies 
into 7 morphology bins: E, SO (including SO and SOa), Sa (Sa 
and Sab), Sb, Sbc, Sc (Sc and Scd), and Sd (13 Sd, 1 Sm and 1 
Irr). 

Fig. 1 shows that these 300 galaxies provide a fair repre¬ 
sentation of the CALIFA survey as a whole. The left panel 
shows scaled histograms of the Hubble type in the mother sam¬ 
ple (empty bars) and in our sample (filled bars). The number of 
objects in each morphology bin for our sample is indicated at 
the top, with a brown to blue color palette that represents Hub¬ 
ble types from ellipticals to late spirals. This same color scheme 
is used throughout this paper. The similarity of the distributions 
reflects the random sampling strategy of CALIFA, with targets 
being picked from the mother sample on the basis of visibility 
criteria alone. The right panel in Fig. 1 shows the u-r versus M r 
CMD, with grey points representing the mother sample and col¬ 
ored points the 300 galaxies. As for the Hubble type distribution. 


a simple visual inspection shows that our subsample is represen¬ 
tative of the full CALIFA sample in terms of CMD coverage. 

2.2. Observations and data reduction 

The observations were carried out with the Potsdam Multi- 
Aperture Spectrometer (Roth et al. 2005, PMAS,) in the PPaK 
mode (Verheijen et al. 2004) at the 3.5m telescope of Calar Alto 
observatory. PPaK contains 382 fibers of 2.7" diameter each, 
and a 74" X 64" Field of View (FoV Kelz et al. 2006). Each 
galaxy is observed with two spectral settings, V500 and V1200, 
with spectral resolutions ~ 6 A (FWHM) and 2.3 A, respectively. 
The V500 grating covers from 3745 to 7300 A, while the V1200 
covers 3650-4840 A. Detailed descriptions of the observational 
strategy and of the data can be found in Sanchez et al. (2012), 
and Husemann et al. (2013). 

The datacubes analyzed here have been calibrated with ver¬ 
sion 1.5 of the reduction pipeline. The main issues addressed 
by this new version are: (i) correction of the sensitivity curve 
for the V500 grating; (ii) new registering method to determine, 
for each galaxy, the relative positioning of the 3 pointings of 
the dithering pattern, and absolute WCS registration; (Hi) a new 
cube interpolation method. CALIFA pipeline vl.5 improves the 
flux calibration to an accuracy of 2-3% and is the current official 
data release. A detailed account of this new pipeline is presented 
in the Data Realease 2 article (Garcfa-Benito et al. 2015). 

In order to reduce the effects of vignetting on the data, we 
combine the observations in the V1200 and V500 setups. The 
combined datacubes were processed as described in Cid Fernan¬ 
des et al. (2013). Our analysis requires that spectra have sig¬ 
nal to noise ratio S/N > 20 in a 90 A window centered at 5635 
A (rest-frame). When individual spaxels do not meet this S/N 
threshold, they are coadded into Voronoi zones (Cappellari & 
Copin 2003). Further pre-processing steps include spatial mask¬ 
ing of foreground/background sources, rest-framing and spectral 
resampling. The resulting 253418 spectra were then processed 
through starlight and pycAsso (the Python CALIFA starlight 
Synthesis Organizer), producing the stellar population properties 
discussed here as described in detail in the next section. 


3. Stellar population analysis: Differences with 
respect to previous work 

Our method to extract stellar population properties from dat¬ 
acubes has been explained and applied to CALIFA in Perez et al. 
(2013), Cid Fernandes et al. (2013, 2014), and Gonzalez Del¬ 
gado et al. (2014c,a). In short, we analyse the data with the 
starlight code (Cid Fernandes et al. 2005), which fits an ob¬ 
served spectrum (0,\) in terms of a model (Mj) built by a non- 
parametric linear combination of V* Simple Stellar Populations 
(SSPs) from a base spanning different ages ( t ) and metallici- 
ties (Z). Dust effects are modeled as a foreground screen with 
a Cardelli et al. (1989) reddening law with Ry = 3.1. Windows 
around the main optical emission lines and the Nal D absorption 
doublet (because of its interstellar component) are masked in all 
spectral fits 2 . Bad pixels (identified by the reduction pipeline) 


2 To test the effect of this process in the estimation of ages, we have 
compared the results for 60 galaxies in common with Sanchez-Blazquez 
et al. (2014). This work uses Steckmap (Ocvirk et al. 2006) and the H J3 
line (previously corrected for emission). Statistically, we find that there 
is no difference in ages (mean = -0.04, std = 0.15 dex) if the same SSP 
models are used in the two methods. 
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Fig. 1 . Left: Comparison of the distribution of Hubble types in the CALIFA mother sample (empty bars) and the galaxies analyzed here (filled 
bars). The number of galaxies in our sample are labeled in colors. The histograms are normalized to form a probability density, i.e., each bar scales 
with the ratio of the number of galaxies in each bin and the total number of galaxies, such that the two distributions are directly comparable. Right: 
Color-magnitude diagram. Mother sample galaxies are plotted in grey, while the 300 galaxies analyzed in this work are marked as colored points. 


are also masked. Results for each spectrum are then packed and 
organized with the pycAsso pipeline. 

This working scheme is preserved here, but with three new 
developments: 

1. The datacubes used in this paper come from the version 1.5 
of the reduction pipeline (Garcla-Benito et al. 2015). 

2. Larger and more complete SSP bases are employed. 

3. A somewhat different definition of mean stellar metallicity is 
adopted (see Gonzalez Delgado et al. 2014a). 

This section describes the novelties related with the stellar 
population synthesis. Improvements resulting from the new re¬ 
duction pipeline are described in Appendix A. 

3.1. SSP spectral bases 

SSP models are a central ingredient in our analysis, linking the 
results of the spectral decomposition to physical properties of 
the stellar populations. Our previous applications of starlight to 
CALIFA explored spectral bases built from three sets of SSP 
models, labeled as GM, CB and BC in Cid Fernandes et al. 
(2014). The first two are again used in this study, but extended 
to a wider range of metallicities, producing what we will denote 
as bases GMe and CBe. 

Base GMe is a combination of the SSP spectra provided by 
Vazdekis et al. (2010) for populations older than t = 63 Myr and 
the Gonzalez Delgado et al. (2005) models for younger ages. The 
evolutionary tracks are those of Girardi et al. (2000), except for 
the youngest ages (1 and 3 Myr), which are based on the Geneva 
tracks (Schaller et al. 1992; Schaerer et al. 1993; Charbonnel 
et al. 1993). The IMF is Salpeter. In our previous studies of the 
first 100 CALIFA galaxies we defined base GM as a regular (f, Z) 
grid of these models, with 39 ages spanning t = 0.001-14 Gyr 
and four metallicities from 0.2 to 1.5 Z Q . We now extend the 
Z range to use of all seven metallicites provided by Vazdekis 
et al. (2010) models: logZ/Z Q = -2.3, -1.7, -1.3, -0.7, -0.4, 0, 
and +0.22. Because these models lack ages below 63 Myr, these 
young ages are only covered by the four largest metallicities, 
such that our extended GM base is no longer regular in t and Z. 
Base GMe contains V* = 235 elements. 


Base CBe is built from an update of the Bruzual & Char¬ 
iot (2003) models (Chariot & Bruzual 2007, private communi¬ 
cation), replacing STELIB (Le Borgne et al. 2003) by a com¬ 
bination of the MILES (Sanchez-Blazquez et al. 2006; Falcon- 
Barroso et al. 2011) and granada (Martins et al. 2005) spectral 
libraries (the same ones used in base GMe). The evolutionary 
tracks are those collectively known as Padova 1994 (Alongi et al. 
1993; Bressan et al. 1993; Fagotto et al. 1994a,b; Girardi et al. 
1996). The IMF is that of Chabrier (2003). Whereas in previous 
works we limited the Z range to > 0.2 solar, we now extend this 
base to six metalicities: log Z/Z 0 = -2.3, -1.7, -0.7, -0.4, 0, 
and +0.4. Base CBe contains A* = 246 elements (41 ages from 
0.001 to 14 Gyr and the 6 metallicities above). 

The main similarities and differences between bases GMe 
and CBe are the same as between the original GM and CB bases, 
thoroughly discussed in Cid Fernandes et al. (2014). Throughout 
the main body of this paper we focus on results obtained with 
base GMe, but we anticipate that our overall qualitative findings 
remain valid for base CBe. The role of base CBe in this paper 
is to allow a rough assessment of the uncertainties associated to 
model choice. 

A minor technical difference with respect to our previous 
analysis is that we now smooth the spectral bases to 6 A FWHM 
effective resolution prior to the fits. This is because the kinemat- 
ical filter implemented in starlight operates in velocity-space, 
whereas both CALIFA and the SSP model spectra have a con¬ 
stant spectral resolution in /(-space, so that effects of the in¬ 
strumental broadening can only be mimicked approximately by 
starlight. We have verified that this modification does not affect 
the stellar population properties used in this paper. 

Appendix B presents some comparisons of the results ob¬ 
tained with these two bases. Experiments were also performed 
with bases which extend the age range to 18 Gyr, and configur¬ 
ing starlight to allow negative values of Ay. These tests are also 
discussed in Appendix B, which adds to the collection of “sanity 
checks" on the results of our analysis. 

4. Galaxy mass, metric, and stellar metallicity 

This section addresses three relatively unrelated aspects, which 
are all important to better understand the results presented in 


Article number, page 4 of 42 


















IAA et al.: Metallicity, age, and stellar mass density across the Hubble sequence 


Table 1. Number of galaxies for each Hubble type and M* interval 
(GMe) 


log M* (M g ) bin 

E 

SO 

Sa 

Sb 

She 

Sc 

Sd 

<9.1 

- 

- 

- 

- 

- 

2 

2 

9.1-9.6 

- 

- 

- 

- 

- 

9 

8 

9.6-10.1 

- 

- 

- 

- 

2 

10 

5 

10.1-10.6 

- 

- 

7 

11 

16 

21 

- 

10.6-10.9 

3 

8 

9 

14 

21 

4 

- 

10.9-11.2 

8 

14 

22 

17 

16 

3 

- 

11.2-11.5 

17 

8 

13 

10 

3 

1 

- 

11.5-11.8 

12 

2 

- 

1 

- 


- 

>11.8 

1 

- 

- 

- 

- 


- 

total 

40 

32 

51 

53 

58 

50 

15 



the next section, where we examine how the spatial distribu¬ 
tion of stellar population properties relates to a galaxy’s stellar 
mass and morphology. First, §4.1 reviews the relation between 
stellar mass and morphological type for our sample. This strong 
relation is imprinted on virtually all results discussed in §5. Sec¬ 
ondly, §4.2 compares our measurements of the Half Light (HLR) 
and Half Mass Radii (HMR). As discussed by Gonzalez Del¬ 
gado et al. (2014c), these two natural metrics for distances are 
not identical due to the inside-out growth of galaxies. Here we 
inspect how the HMR/HLR ratio varies as a function of Hubble 
type and stellar mass in our sample. Finally, §4.3 presents our 
definition of mean stellar metallicity. Gonzalez Delgado et al. 
(2014c) showed that stellar mass surface densities, mean ages, 
and extinction values defined from the integrated spectrum, from 
galaxy-wide spatial averages, and measured at R = 1 HLR all 
agree very well with each other. Here we extend this test to stel¬ 
lar metallicities. Throughout this section, results for the two SSP 
models discussed in §3.1 are presented. 


4.1. Stellar masses 

To obtain the total stellar mass of a galaxy we add the mass in 
each zone, thus taking into account spatial variations of the stel¬ 
lar extinction and M/L ratio. Masked spaxels (e.g., foreground 
stars) are accounted for using the /i* radial profile as explained 
in Gonzalez Delgado et al. (2014c). 

Fig. 2 shows the distribution of M* as a function of Hub¬ 
ble type. Table 1 shows the distribution of galaxies by Hubble 
type in several bins of M*. The masses range from 7 x 10 8 to 
7 x 10 n M Q for fits with GMe (Salpeter IMF). CBe-based masses 
(Chabrier IMF) are on average smaller by a factor 1.84. As for 
the general galaxy population, mass is well correlated with Hub¬ 
ble type, decreasing from early to late types. High bulge-to-disk 
ratios (E, SO, Sa) are the most massive ones (> 10 n M o ), while 
galaxies with small bulges (Sc-Sd) have M* < 1O 1O M 0 . The av¬ 
erage log M*(M q ) is 11.4, 11.1, 11.0, 10.9, 10.7, 10.1, and 9.5 
for E, SO, Sa, Sb, She, Sc, and Sd, respectively. The dispersion 
is typicaly 0.3 dex, except for Sc galaxies, that have a dispersion 
of ~ 0.5 dex. 

Because CALIFA is not complete for M, > -19.5, this dis¬ 
tribution in mass is not completely representative of the local 
Universe. In particular, it is important to remember that dwarf el¬ 
lipticals are not included, so M* or any other property discussed 
here for E’s are restricted to massive ellipticals. 


Fig. 2. Distribution of the stellar masses obtained from the spatially 
resolved spectral fits of each galaxy for each Hubble type (grey small 
points). The colored dots (stars) are the mean galaxy stellar mass in each 
Hubble type obtained with the GMe ( CBe) SSP models. The bars show 
the dispersion in mass. 

4.2. The HMR/HLR 

As explained in Cid Fernandes et al. (2013), we define the HLR 
as the semi-major axis length of the elliptical aperture that con¬ 
tains half of the total light of the galaxy at the rest-frame wave¬ 
length 5635 A. Similarly, the HMR is derived from the 2D dis¬ 
tribution of the stellar mass, as the elliptical aperture at which 
the mass curve of growth reaches 50% of its asymptote. The ra¬ 
tio between the HMR and the HLR (ci^Jci^ Q ) reflects the spatial 
variation of the star formation history in a galaxy. This ratio is 
lower than 1 in almost all cases (Gonzalez Delgado et al. 2014c), 
a signpost of the inside-out growth found by Perez et al. (2013). 

Fig. 3 shows the relation between and Hubble type 

(left panel), and galaxy stellar mass (right panel). These plots 
confirm our earlier finding that galaxies are more compact in 
mass than in light. If the gradient in stellar extinction is taken 
into account, the average = 0.82(0.80) + 0.10(0.13) 

for base GMe (CBe). Fig. 3 shows that the ratio decreases from 
late to early type spirals; while lenticulars and ellipticals have 
similar a^ Q /a^ Q . 

These results are also in agreement with our previous result 
that a^Ja^Q shows a dual dependence with galaxy stellar mass: It 
decreases with increasing mass for disk galaxies but it is almost 
constant in spheroidal galaxies, as confirmed in the right panel 
of Fig. 3. Sb-Sbc galaxies are the ones with the lowest af Q /a^ iy 

4.3. Stellar metallicity 

Metallicity is one of the most difficult stellar population proper¬ 
ties to estimate. Reasons for this include: (i) the coarse metal¬ 
licity grid of the SSP bases; (ii) the limitation of the stellar li¬ 
braries to the solar neighborhood; and (Hi) inherent degeneracies 
like the dependence of the continuum shape on extinction, age, 
and metallicity, whose effects are hard to disentangle. Notwith¬ 
standing these difficulties, meaningful estimates of Z* can be 
extracted from observed spectra, particularly by means of full 
spectral synthesis methods (Sanchez-Blazquez et al. (2011)). 

STARLiGHT-based estimates of Z* for the same CALIFA sam¬ 
ple used in this paper were previously used by Gonzalez Del- 
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Fig. 3. Left: The ratio between half mass and half light radius (a^>/ aj 0 ) with the Hubble type (left). Big colored dots represent the averaged a^,/ aj 0 
in each Hubble type bin, and the lines the dispersion. Stars and big circles show the results obtained with the GMe and CBe bases, respectively. 
Right: a^,/ af 0 as a function of the galaxy stellar mass. The black circles show the averaged correlation independently of the morphological type. 
Large circles represent the averaged relation in mass intervals of 0.25 dex for each color-coded morphological type. 


gado et al. (2014a) to study global and local relations of Z* with 
the stellar mass and stellar mass surface density. We have shown 
there that: (i) our sample follows a well defined stellar mass- 
metallicity relation (MZR), (ii) this relation is steeper than the 
one obtained from O/H measurements in HII regions, but that 
considering only young stellar populations the two MZR’s are 
similar, and (iii) Z* is strongly related to /r* in galaxy disks and 
to M* in spheroids. All these results lend confidence to our Z* 
estimates. 

Here we review our definition of the mean stellar metallic- 
ity, and test whether its value at 1 HLR matches the galaxy wide 
average value as well as the one obtained from the spatially col¬ 
lapsed data cube. 

4.3.1. Mean stellar metallicity 

The main properties analyzed in this paper are the stellar 
mass surface density (//*), stellar extinction (Ay), mean age 
((log age)i ), and metallicity of the stellar population, whose spa¬ 
tial distributions are studied as a function of Hubble type and to¬ 
tal stellar mass (M*). These properties were defined in previous 
articles in this series. For instance, the mean light weighted log 
stellar age is defined as 

(log age) L = 'Yj x tz x log t (1) 

t,z 

(eq. 9 of Cid Fernandes et al. 2013), where x,/ is the fraction of 
flux at the normalization wavelength (5635 A) attributed to the 
base element with age t and metallicity Z. The mass weighted 
version of this index, (log age)M , is obtained replacing x,/ by its 
corresponding mass fraction m t/ . 

While Cid Fernandes et al. (2013) average the base metallic- 
ities linearly (their eq. 10), in this paper, as in Gonzalez Delgado 
et al. (2014a), we employ a logarithmic average: 

(log Z*) M = 'Yj m tz x log Z (2) 

t,Z 

for the mass weighted mean log Z* and 


(log Z*)l = ^ x,z x log Z (3) 

t,z 

for the luminosity weighted mean logZ*. The motivation to use 
this definition is that the extended SSP bases used in this study 
span a much wider dynamical range in Z* (nearly three orders 
of magnitude, compared to barely one in our previous papers), 
which is better handled with a geometric mean (implicit in the 
use of the logarithm). This is the same reasoning behind the use 
of (log age) instead of log {age). 

To some degree, the definition of mean Z is largely a matter 
of taste (albeit one with mathematical consequences because of 
the inequality of the arithmetic and geometric means, (log Z) < 
log(Z)), so much so that one finds both types of averaging in the 
literature. For instance, in Gallazzi et al. (2005) metallicities are 
averaged logarithmically, whereas Asari et al. (2007) work with 
arithmetic averages. 

As shown in Gonzalez Delgado et al. (2014a) (see also Fig. 
4), our metallicities span about 1 dex for galaxy masses ranging 
from 10 9 to 10 12 M g , with an MZR which matches well the stellar 
metallicities of both Milky Way and LMC-like galaxies. 

4.3.2. Galaxy averaged stellar metallicity 

Gonzalez Delgado et al. (2014c) obtained the important result 
that galaxy-averaged stellar ages, mass surface density, and ex¬ 
tinction are well matched by the corresponding values of these 
properties at R = 1 HLR and also with the values obtained from 
the analysis of the integrated spectrum (i.e, the one obtained by 
collapsing the datacube to a single spectrum). The general pat¬ 
tern therefore is that galaxy averaged properties match both the 
values at 1 HLR and those obtained from integrated spectra. Do 
our stellar metallicities comply with this rule? 

To answer this question we first define the galaxy-wide av¬ 
erage stellar metallicity following eq. 2 in Gonzalez Delgado 
et al. (2014a). which gives the mass weighted mean value of 
(log Z* >;cy > M as 

„_ 7 x galaxy _ Sat 4^*,xy(log Z+) M ,xy 

U°g Z *>M - V AA 
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Fig. 4. Upper panels: Comparison of the galaxy-wide average stellar 
metallicity (weighted in mass) derived from the spatially resolved spec¬ 
tral analysis ((logZ*)^ ,ux> ) and the integrated metallicity derived from 
fitting the total (integrated) galaxy spectrum ((log Z i ,) , ^ esm ‘ ed ). Lower 
panel: Comparison of (log Z-i,) s ^ laxy with the value measured R = 1 
HLR ((log Z*)^ LR ). Left and right panels show results obtained with 
base GMe and CBe SSPs, respectively. All panels include 300 galaxies. 
The difference between the y-axis and x-axis is labeled in each panel as 
A, and the dispersion as cr. 


where M * iXy is the stellar mass in spaxel xy. 

Fig. 4 compares our results for (logZ*)^ /fl ' v ' with the mass 
weighted mean (log Z*) ;M values obtained at R - I HLR 
((log Z*)^ LR , bottom panels) and those derived from the inte¬ 
grated spectrum ((log Z- k )"^ egra,ed , top panels), analyzed in the 
exact same way as the individual zone spectra. Results are shown 
for both base GMe (left panels) and CBe (right). 

The agreement is remarkable. The galaxy averaged metallic¬ 
ity and the one at 1 HLR are the same to within a dispersion of 
0.1 dex. The integrated metallicity also matches the galaxy aver¬ 
aged value, with only slightly larger dispersions. The largest de¬ 
viations occur for low metallicity systems. Similar conclusions 
are reached if the comparison in Fig. 4 is done using the light 
weighted version of eq. (4), 


<logZ.)—” = (5) 

2axy B-k,xy 

where L* vv is the luminosity (corrected by stellar extinction) in 
each spaxel evaluated at a reference wavelength (5635 A in our 
case). 

The stellar metallicities behave as expected, in the sense 
that, like other properties, their galaxy-wide averages match the 
values at R — 1 HLR, and also the values derived from in¬ 
tegrated spatially unresolved spectroscopy (Gonzalez Delgado 
et al. 2014c). We thus conclude that galaxy-wide spatially av¬ 
eraged stellar population properties (stellar mass, mass surface 
density, age, metallicity, and extinction) match those obtained 
from the integrated spectrum, and that these spatially averaged 
properties match those at R= 1 HLR, proving that effective radii 
are really effective (Gonzalez Delgado et al. 2014b). 


5. Spatially resolved stellar population properties 
as a function of morphology and mass 

This section presents a series of results derived from our spa¬ 
tially resolved spectral synthesis analysis of CALIFA galaxies. 
We focus on the following four stellar populations properties: 
mass surface density (p t , §5.1), mean ages ((log age) L , §5.2), 
metallicities ((logZ*)^, §5.3), and extinction (Ay, §5.4). Each 
of these properties is studied by means of (i) 2D maps of the 
individual galaxies, (ii) radial profiles, and (Hi) radial gradients. 
Throughout the section, the emphasis is on evaluating and com¬ 
paring the roles of morphology and total stellar mass in shaping 
the observed behavior of these four properties. 

Before discussing the results, we briefly explain how these 
quantities are obtained and how they are presented. 

2D maps in the CMD: Using pycAsso we obtain, for each galaxy, 
2D maps of each of the four properties. The results for all the 
galaxies are presented in the framework of the color-magnitude 
diagram, where each map is placed at the galaxy’s coordinates 
in the m - r vs. M r CMD. Because absolute magnitude is re¬ 
lated to M* and redder galaxies are (usually) older and more 
metal rich, these plots show the correlations M*-age, 

and M*-metallicity in a 2D fashion. Because in our sample the 
galaxy Hubble type is correlated with color and luminosity, these 
plots not only show how the galaxy averaged properties and their 
radial structure change with the galaxy stellar mass, but also with 
the morphological type. These maps are shown in the Appendix 
C (Figs. C.1-C.4). 

Radial profiles: Each 2D map is azimuthally averaged in order 
to study the radial variations of each of the four stellar popula¬ 
tion properties. Elliptical apertures 0.1 HLR in width are used to 
extract the radial profiles, with ellipticity and position angle ob¬ 
tained from the moments of the 5635 A flux image. Expressing 
radial distances in units of HLR allows the profiles of individual 
galaxies to be compared on a common metric, as well as aver¬ 
aging (“stacking”) them as a function of Hubble type or stellar 
mass. Radial profiles expressed in units of the HMR were also 
analyzed and lead to similar shapes, so all profiles presented be¬ 
low use HLR as the unit for radius. 

Radial gradients: Inner and outer gradients are defined as differ¬ 
ences between the values at R — 1 and 0 (V,„), and R = 2 and 1 
(V ONt ), respectively. For instance, 

V,„ logp* = log/.<*(! HLR) - log/r*(0) (6) 


Vou, logp* = logp*(2 HLR) - logM1 HLR) (7) 


for logp*, and similarly for (log age)i, (log Zfi) M and A v . De¬ 
fined in this way, the gradients have units of dex/HLR (mag/HLR 
for vAy). Since the stellar population properties of galaxies at 1 
HLR represent very well the galaxy-wide average, V,„ (V OH ,) ef¬ 
fectively measures how the bulge (disk) properties change with 
respect to those of the galaxy as a whole. 3 

3 Based on the exponential fit analysis developed by (Sanchez et al. 
2013) and (Sanchez-Blazquez et al. 2014), we conclude that the regions 
between 1 and 2 HLR are dominated by the disk component; thus, V ouf 
measures the disk gradient. However, V,„ is not measuring the bulge 
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Unless otherwise noted, all results reported below are for 
the GMe base, although the whole analysis was carried out with 
properties derived with both sets of SSP models discussed in 

3.1. Differences between GMe and CBe SSPs go in the follow¬ 
ing way: (a) The stellar mass surface density is lower with CBe 
than with GMe by 0.27 dex on average, mostly due to the differ¬ 
ent IMFs (Salpeter in GMe versus Chabrier in CBe). (b) Varia¬ 
tions in stellar extinction are negligible, (c) CBe yields somewhat 
younger ages and higher metallicities than GMe, by an average 
of 0.14 dex in (log age)L and 0.12 dex in (logZ*)^. These shifts 
reflect the age-metallicity degeneracy, and are mainly a conse¬ 
quence of the different sets of metallicities available in these 
bases. However, radial gradients are not affected by this degen¬ 
eracy. A detailed comparison of properties derived with the two 
bases is given in Appendix B. 

5.1. Stellar mass surface density 

2D maps of the stellar mass surface density for the 300 individual 
galaxies of our sample are presented in the Appendix C (Fig. 
C.l). Here we discuss the radial structure of log/r* as a function 
of Hubble type and M*. 

5.1.1. yu*-morphology and /r*-mass relations 

Fig. 5 shows how //* measured at 1 HLR changes with Hubble 
type (left panel), and with the galaxy stellar mass (right). Recall 
from Gonzalez Delgado et al. (2014c) that properties measured 
at 1 HLR match very well the corresponding galaxy-wide av¬ 
erage value, so these plots ultimately show how the global /./„ 
depends on the morphology and on M*. 

The plot shows (log/j+ LR ) increasing from late spirals to 
spheroids, with average and dispersion values of 3.1 ± 0.2, 
3.10 ± 0.18, 3.05 ± 0.25, 2.70 + 0.17, 2.65 + 0.24, 2.40 ± 0.28, 
2.04 + 0.27, for E, SO, Sa, Sb, She, Sc and Sd, respectively. Note 
that E and SO are remarkably similar. 

Surface densities also increase with M*, as seen in the 
right panel of Fig. 5. The overall relation is relatively 

smooth, with no evidence of an abrupt change of behavior as 
that discussed by Kauffmann et al. (2003) for SDSS galaxies. 
Fig. 5, however, reveals that morphology is also behind the dis¬ 
persion in the relation. The black line shows the relation 

for the full sample, obtained by averaging log /.(* in 0.4 dex-wide 
bins in mass, while the big circles break this general relation into 
different (color-coded) morphological types for the same mass 
bins. Despite the reduced statistics, it is evident that: (a) for the 
same stellar mass, early type galaxies are denser than late type 
ones, and (b) Sa and earlier type galaxies exhibit a much flat¬ 
ter /r*-M* relation than later types. The overall impression from 
these results is that morphology, and not only stellar mass, plays 
a fundamental role in defining stellar surface densities, and it is 
responsible for the change of slope in the SDSS relation. 

5.1.2. Radial Profiles 

Azimuthally averaged radial profiles of log/r* are shown in Fig. 
6. Results are stacked by Hubble type (left panel) and mass 
(right). In the left panel galaxies are grouped in our seven mor¬ 
phological classes. The typical dispersion within these bins is 

gradient. The reason is that the effective radius (R f ) of the spheroidal 
component can be smaller than 1 HLR, and it shows a dependence with 
the morphological type. Thus, R f ~ 1 HLR for E, but is significantly 
smaller in late type spirals. 


illustrated by the error bar, which shows the standard deviation 
in log /r* (R - 1 HLR) for galaxies of the Sa class. 

A clear trend with Hubble type is seen: The p+(R) profiles 
scale with Hubble type from late to early spirals, and this mod¬ 
ulation with morphology is preserved at any given distance. E 
and SO have remarkably similar profiles, with core and extended 
envelope equally dense at any given distance, suggesting that the 
disk of SO galaxies and the extended envelope of ellipticals have 
grown their mass by similar processes. 

The right panel of Fig. 6 shows the radial profiles grouped in 
seven bins of stellar mass spanning the log M*(M 0 ) = 9.1-11.8 
range. These also show that the average of log /r*(R) is mod¬ 
ulated by M*. However, this /r*(R)-M* modulation breaks for 
early type galaxies (concentration index C (Vyo/rso) >2.8; see 
also Fig. 13 in Gonzalez Delgado et al. (2014c)), that in our sam¬ 
ple are populated mainly by E and SO, and some Sa. On the other 
hand these early types are all massive, with M* > 10 11 M G . 

5.1.3. Radial gradients 

Inner (0-1 HLR) and outer (1-2 HLR) gradients in log//*, as de¬ 
fined by equations (6) and (7), are plotted as a function of mor¬ 
phology and stellar mass in Fig. 7. V,„ log /r* values (correspond¬ 
ing to the core region) are plotted in grey-red, while V out log /r* 
(which trace the disks of spirals and SO and the extended en¬ 
velope of ellipticals) are plotted in grey-blue. Circles and stars 
show results for bases GMe and CBe respectively, illustrating 
that even though these bases yield different absolute values of 
/.(* the resulting gradients are nearly identical. 

A clear correlation exists between V,„log/r* and Hubble 
type. The gradient in the inner HLR increases (in absolute val¬ 
ues) significantly from late to early spirals, converging to a con¬ 
stant value for E and SO. This relation reflects the variation of 
the bulge to disk ratio in spirals, and the dominance of the bulge 
component in spheroids (SO and E). The outer gradient is weaker 
(smaller in absolute value) than the inner one, as expected if a 
disk component dominates the mass outwards of 1 HLR. 

The right panel of Fig. 7 shows the relation between the in¬ 
ner gradient and the stellar mass. There is a clear increase (in 
absolute values) of V,„ log/r* with M*, with the more massive 
galaxies having a steeper increase of the central density. The 
dispersion with respect to the average values (black cross) within 
M*-bins is significant. To check the effect of morphology on this 
dispersion we have averaged V,„ log /r* in mass intervals for each 
Hubble type and plotted the resulting averages (large colored cir¬ 
cles). The general trend that emerges is that, for galaxies of the 
same mass, early type galaxies tend to be overall centrally denser 
than later types, in agreement with Fig. 5; although, there are a 
few intervals of stellar mass (e.g. log M* = 11.4 M 0 ), in which 
the variations in V,„ log /r* with Hubble type are not significant. 

It is also worth mentioning that V,„ login Sa and Sb is 
very close to that in SO and E, and in this sense it would be easy 
to fade early type spirals into SO’s. 

5.2. Ages of the stellar populations 

2D maps of the luminosity weighted mean log stellar ages (eq. 
1) for the 300 galaxies are presented in Fig. C.2. Here we discuss 
the radial structure of (log age)i and its relation to Hubble type 
and M*. The presentation follows the same script used in the 
presentation of //^-related results in §5.1. 
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Fig. 5. Left panel: stellar mass surface density measured at 1 HLR as a function of Hubble type. Small dots represent log/j* for each galaxy; the 
colored circles are the average logyu* for each Hubble type, and the error bars are the dispersion in log/j* for each morphological type. Right 
panel: log/r*-log M* relation. Individual galaxies are represented by small dots colored by their morphological type. The black line is the average 
log/r* in galaxy stellar mass bins of 0.4 dex. Large colored circles are the average logp* in each bin of mass for each Hubble type. 




Fig. 6. (left) Radial profiles (in units of HLR) of the stellar mass surface density obtained with base GMe. The results are stacked in seven 
morphology bins. The error bar in the panel indicate the dispersion at one HLR distance in the galaxies of the Sa bin. It is similar for other 
Hubble types and radial distances, (right) Radial profiles stacked in seven bins of galaxy stellar mass, log M*(M G ): 9.1-9.6, 9.6-10.1, 10.1-10.6, 
10.6-10.9, 10.9-11.2, 11.2-11.5, 11.5-11.8. 


5.2.1. Age-morphology and age-mass relations 

Fig. 8 shows how the mean age of the stellar populations 
at 1 HLR changes along the Hubble sequence (left panel), 
and with the galaxy stellar mass (right). Similarly to log/r* LS , 
(log age)i HLR represents well the galaxy-wide averaged stel¬ 
lar population age ((log age)i, galaxy , Gonzalez Delgado et al. 
(2014c)). 

Clearly, (log age)i HLR scales with Hubble type, increasing 
steadily from Sd to Sa. SO and ellipticals have stellar populations 
of similar mean age, and older than spirals. The average and dis¬ 
persion values of (log age)L HLR (yr) are 8.62 ± 0.22, 8.89 ± 0.22, 
9.07 + 0.19, 9.33 + 0.21,9.55 + 0.19, 9.71 ±0.11, and 9.74 + 0.11, 
for Sd, Sc, She, Sb, Sa, SO and E, respectively. 

Mean ages also increase with the galaxy mass (right panel of 
Fig. 8), a “downsizing" behavior that has been confirmed with 
widely different samples and methods. For instance, our age- 
mass relation is similar to that derived for SDSS galaxies by 
Gallazzi et al. (2005) (their figure 8). They found that there is 


a transition at M* ~ 3 x 1O 1O M 0 4 , below which galaxies are 
typically young and above which they are old. This is the same 
mass at which Kauffmann et al. (2003) find the relation 

to flatten. 

Unlike in these SDSS-based works, we do not see sharp tran¬ 
sitions as a function of M* in neither //* nor (log age) L , although 
differences in sample selection and statistics prevent a proper 
comparison. We do, however, note a common behavior in the 
right panels of Figs. 5 and 8, in the sense that the dispersion 
above ~ 1O 1O M 0 is strongly related to morphology. 

Like in Fig. 5 (right panel), the black line in the right panel 
of Fig. 8 shows the age-mass relation for the whole sample, ob¬ 
tained by averaging (log age)i HLR values in M* bins. Small dots 
show individual galaxies, while the large colored circles rep¬ 
resent the mass-binned average (log age)L HLR for each Hubble 
type. As with the /r*-M* relation, breaking the age-mass rela¬ 
tion into morphological types reveals clean trends. In this case, 

4 Equivalent to ~ 5.5 x 10 10 M o for our IMF. 
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Fig. 7. (left) Correlation between the inner (grey-red) and outer (grey-blue) gradient of log /j* and the morphological type. The results are shown 
for the GMe (stars) and CBe (circles) SSP models. The inner gradient is calculated between the galaxy nucleus and 1 HLR, and the outer gradient 
between 1 HLR and 2 HLR. (right) Correlation between the inner gradient of log /r* and the galaxy stellar mass. Small dots represent the results 
for each galaxy, and black crosses the average for each 0.3 dex mass bin. Large circles represent the averaged inner gradient in mass intervals of 
0.3 dex for each color-coded morphological type. Black crosses show the average correlation between the inner gradient of log yu* and galaxy mass 
independently of the morphological type. 


we see that, for a fixed M*, earlier type galaxies are older than 
later types. The corollary is that mass is not the sole property 
controlling the SFH of a galaxy. In fact, given the ~ flat age- 
mass relations for Sa, SO and E, morphology seems to a more 
relevant factor, at least in these cases. 

5.2.2. Radial profiles 

Fig. 9 shows the age radial profiles obtained by stacking galax¬ 
ies as a function of Hubble type and mass. The {\og age) l(R) 
profiles scale with Hubble type, but by different amounts at the 
center than at 1 HLR. At any radial distance, however, the early 
type galaxies are older than later type ones. E and SO are again 
very similar at all radii. This suggests that E and SO have simi¬ 
lar histories not only on average, but also spatially resolved, at 
least in the inner 2 HLR. Negative age gradients are detected in 
all galaxies (except perhaps in Sd, whose ages profiles are flat¬ 
ter than in the other spirals 5 ). These negative gradients reflect 
the inside-out growth of galaxies. Furthermore, the decrease of 
(log age)i with R indicates that quenching happens earlier at the 
galaxy center; and also earlier in early type galaxies (spheroids 
and Sa) than in later type spirals (Sbc-Sc). 

The radial profiles also show a clear trend with M* (Fig. 9, 
right), with the more massive galaxies being older everywhere, 
hence preserving the downsizing pattern at all radial distances. 
Comparing the left and right panels in Fig. 9, one sees that 
grouping galaxies by their stellar mass leads to a reduced ver¬ 
tical stretch in their (log age)i(R) profiles than when the aver¬ 
aging is done by morphology. But the profiles expand similar 
vertical scale if galaxies earlier than Sd and more massive than 
10 9 6 M 0 are considered; indicating that the effect of morphol¬ 
ogy and stellar mass are not easily disentangled here. However, 
in §6.3, Fig. 20 shows that the dispersion in the (log age)r(R) 
profiles between galaxies of the same M* and different Hubble 


5 The small drop of (log age)i toward the center of Sd galaxies is 
caused by a couple of galaxies with young nuclear regions. Given that 
this group is the least populated in our analysis (only 15 galaxies), better 
statistics is needed to evaluate the reality of this feature. 


type is significant, and larger than between the (log age)i(R) pro¬ 
files of galaxies of different M* but the same Hubble type. These 
results in agreement with Fig. 8 indicate that the age profiles are 
more related to morphology than to M*. Since (log age)/ (R) is 
essentially a first moment of the spatially resolved SFH, we can 
conclude that the SFH and its radial variation are modulated pri¬ 
marily by the galaxy morphology, with mass playing a secondary 
role. 


5.2.3. Radial gradients 

Gradients in (log age)i, computed as indicated in eqs. (6) and 
(7), are plotted in Fig. 10 against Hubble type (left panel) and 
stellar mass (right). The figure layout is exactly as in Fig. 7. 
Whilst in that plot results obtained with bases GMe and CBe 
(circles and stars in the left panel, respectively) could hardly be 
distinguished, here the results for these two sets of SSPs do not 
overlap so precisely, although the differences in V(log age)/ are 
clearly very small (see §B.2). 

A clear relation exists between V,„(log age)r and morphol¬ 
ogy: The inner age gradient increases from early type galaxies 
to Sb-Sbc spirals, which are the galaxies with the largest varia¬ 
tion between the age of the stellar population at the bulge and 
the disk. Spirals of later type (Sc and Sd) have flatter radial pro¬ 
files than Sb-Sbc. The outer (between 1 and 2 HLR) age gradient 
shows a similar bimodal behavior as V,„ (log age)i, but with a 
smaller amplitude. 

The right panel of Fig. 10 shows the behavior of 
V,„(log age)r with M*. The gradient tends to increase (become 
more negative) from low mass galaxies (which have roughly 
flat profiles) up to about 10 n M G , at which point the trend re¬ 
verses and V,„(log age)i decreases with increasing M*. This is 
best seen following the black crosses, that trace the mass-binned 
mean relation. The dispersion with respect to this relation is sig¬ 
nificant and is related to the morphology, as seen through the 
large colored circles. The tendency is that, at a given mass, SO 
and early type spirals have weaker V,„(log age)i than Sb-Sbc. 
This dependence of age gradients with the Hubble type at a fixed 
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Fig. 8. (log age)L measured at 1 HLR as a function of Hubble type (left) or galaxy stellar mass (right). Symbols and colors are as in Fig. 5. The 
black line is the average (log age) L in galaxy stellar mass bins of 0.4 dex. 




0123 0123 

radial distance (HLR) radial distance (HLR) 

Fig. 9. Radial profiles of (log age)L as a function of Hubble type (left) and in seven bins of galaxy stellar mass (right). These bins are log M*(M 0 ) 
= 9.1-9.6, 9.6-10.1, 10.1-10.6, 10.6-10.9, 10.9-11.2, 11.2—11.5, 11.5-11.8. Symbols and colors are as in Fig. 6. These results are obtained 
with base GMe. 




Fig. 10. (left) As Fig. 7 but for (log age)i. The inner gradient shows a clear dependence with Hubble type, that seems to be stronger than with the 
galaxy mass. Sb-Sbc-Sc galaxies have larger inner gradients with CBe than with GMe , but both sets of models show a similar dependence with 
Hubble type, (right) The inner gradient of (log age)i as a function of galaxy mass. Colors and symbols are as in Fig. 7. 
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log M* (M 0 ) 


Fig. 11. (logZ*) M measured at 1HLR as function of Hubble type (left) 
and galaxy stellar mass (right). Symbols and colors are as in Fig. 5. The 
black line is the average (logZ*)iw obtained in 0.4 de bins of log M* . 


M* indicates again that the spatial variation of the SFH is mainly 
driven by the morphology and not by the stellar mass. 

However, the morphology (understood as the B/D ratio (Gra¬ 
ham & Worley 2008)) can not be the only driver of the spatial 
variation of the SFH along all the Hubble sequence. Fig. 10 
shows that there is not a monotonic relation between the B/D 
ratio and V,„(log age)i , with galaxies with the smaller B/D ra¬ 
tio having the largest variations in (log age)i between the central 
core and the disk. This bimodal behavior seen in Fig. 10 suggests 
that other physical properties are also important in establishing 
the spatial variation of the SFH, which on the other hand is re¬ 
flecting the different bulge formation processes along the Hubble 
sequence. 


5.3. Stellar metallicity 

Fig. C.3 presents the images of the mass weighted mean (loga¬ 
rithmic) stellar metallicity (cf. eq. 2). Here we discuss the radial 
structure of (log Z*),w as a function of Hubble type and M*. 


5.3.1. Metallicity-morphology and mass-metallicity relations 

Fig. 11 shows how the stellar metallicity measured at 1 HLR 
changes with the Hubble type (left panel) and with the galaxy 
stellar mass (right). 

Stellar metallicities grow systematically from late to early 
type galaxies. The statistics within each Hubble type are 
(log Z*> M ™*(Z 0 ) = -0.05 ± 0.13, -0.05 ± 0.33, -0.21 + 0.16, 
-0.10 ± 0.18, -0.05 ± 0.15, +0.06 ± 0.08, and +0.10 + 0.08 for 
Sd, Sc, She, Sb, Sa, SO, and E, respectively. 

Not surprisingly, metallicities also grow with M*, as shown 
in the right panel of Fig. 11. Since we have shown in §4.3.2 that 
the galaxy-wide average stellar metallicity is well represented 
by the metallicity at 1 HLR, this plot is in fact equivalent to the 
global mass-stellar metallicity relation (MZR). We have previ¬ 
ously found that this relation is steeper than the one derived from 
HII regions, which is similar to the flatter stellar MZR obtained 
when we consider only young stars (Gonzalez Delgado et al. 
2014a). As in Fig. 5, the smoothed black curve is obtained by av¬ 
eraging (log Zi,)M HLR in 0.4 dex bins of log M*. The dispersion 
in the MZR is significant, and larger than the dispersion pro¬ 
duced by the galaxy morphology as shown by the distribution of 
large colored circles. These circles are the average {\ogZ+) m hlr 
in each mass bin for each Hubble type, and show the tendency of 
earlier type galaxies to be more metal rich than late type galaxies 
of the same stellar mass. 


5.3.2. Radial profiles 

Fig. 12 shows the results of stacking the radial profiles of 
(log Z*)m as a function of Hubble type and M*. Outwards de¬ 
creasing (logZ*)^ is detected for most morphological classes, 
but flat profiles are found for Sc-Sd galaxies. Intermediate type 
spirals (Sb-Sbc) stand out as the ones with the largest variations 
in stellar metallicity. 

The behavior of the radial variation of the stellar metallicity 
with M* (right panel in Fig. 12) is similar to the behavior with 
morphology. Most galaxies have (log Z*) w that decreases with 
R, except for the two lowest mass bins, which show flat profiles. 
The largest spatial variations are also found in galaxies in the 
intermediate mass bins (10 < log M*(M 0 ) < 11). 

These negative radial gradients of the metallicity are also an 
indicator of the inside-out formation processes in galaxies. The 
inversion of the gradient in late type spirals and in low mass spi¬ 
rals may be an indicator of the secular processes or the outside-in 
formation scenario in these galaxies (Perez et al. 2013). 


5.3.3. Radial gradients 

Fig. 13 clones Figs. 7 and 10 for (logZ*)^ gradients. On the 
left panel, one sees that, as for stellar densities and ages, results 
for bases GMe and CBe are very similar. On average, galaxies 

have (log Z») v/ gradients-0.1 dex/HLR, similar to the value 

obtained from nebular oxygen abundances (Sanchez et al. 2013). 
Outer and inner gradients are not significantly different. Despite 
the large scatter, there is a hint of a bimodal distribution as that 
found for stellar ages, also with intermediate type spirals in a 
pivotal position and late type spirals with the flattest gradients, 
at least in a statistical sense. 

The right panel of Fig. 13 shows V,„(logZ*)M as a function 
of M*. The dispersion is significant, but on average there is a 
tendency to turn flat profiles into negative gradient ones as M* 
increases from 10 9 to 10 1() M Q . The largest gradients are found 
between 10 10 and 1O U M 0 . More massive galaxies tend to have 
weaker stellar metallicity gradients. The dispersion is significant 
throughout this relation. A trend with morphology is seen in the 
sense that, for a given mass, early types are the ones with weaker 
gradients. 


5.4. Stellar extinction 

starlight models the stellar extinction as a foreground screen, 
parametrized by Ay and following the Galactic reddening law. 
Images showing the spatial distribution of Ay for our 300 galax¬ 
ies are presented in Fig. C.4. Here we present Ay related results 
as a function of Hubble type and M*, following the same script 
adopted in the discussion of p*, (log age)i. and (logZ *)m in the 
previous subsections, thus completing the list of stellar popula¬ 
tion properties studied in this work. Unlike masses, ages, and 
metallicities, extinction is more easily affected by inclination ef¬ 
fects, so the results reported below should be interpreted with 
caution. Section 5.5 explores this issue in depth. 


5.4.1. Extinction-morphology and extinction-mass relations 

Fig. 14 shows how the stellar extinction at 1 HLR changes with 
Hubble type (left panel), and with stellar mass (right panel). As 
with other properties, Ay LR represents well the mean extinction 
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Fig. 12. Radial profiles of (logZ*) M as a function of Hubble type (left) and of galaxy stellar mass (right). Mass bins are log M*(M 0 ) = 9.1-9.6, 
9.6-10.1, 10.1-10.6, 10.6-10.9, 10.9-11.2, 11.2—11.5, 11.5-11.8. Symbols and colors are as Fig. 6. These results are obtained with base GMe. 




Fig. 14. A v measured at 1HLR as function of Hubble type (left) and 
galaxy stellar mass (right). Symbols and colors are as in Fig. 5. The 
black line is the average (log Z* ) M obtained in 0.4 dex bins of log M*. 


of the galaxy 6 as well as the Ay value derived from spectral fits 


6 The galaxy average extinction for each galaxy is calculated as the 
mean of all the 20 radial values obtained for each galaxy between the 
center and 2 HLR. 


of the integrated spectra. 7 The left panel in Fig. 14 shows Ay LR 
as a function of morphology. Ellipticals and SOs have almost no 
extinction, with mean Ay LR = 0.01 ± 0.01, and 0.06 + 0.07 mag, 
respectively. Sa, Sb and Sc galaxies have Ay LR around 0.25 mag, 
and somewhat smaller (0.19 ± 0.08 mag) in Sd’s. 

There is no clear behavior of stellar extinction with galaxy 
stellar mass. In general, galaxies with M* < 10 n M o have 
Ay LR = 0.2-0.3 mag. More massive galaxies are less extin¬ 
guished, and for fixed mass early types tend to have smaller 
Ay LR , but the dispersion is large. 

5.4.2. Radial profiles 

Fig. 15 shows Ay(R ) profiles stacked by Hubble type (left panel), 
and mass (right). Spirals have A v ~ 0.2 mag in the disk and up to 
0.6 mag in the center. Their Ay profiles are similar for all Hubble 
types except for Sd’s, where Ay does not change much from disk 
to center. Ellipticals and SO’s also show negative Ay gradients, 
although at distances larger than 1 HLR they are almost dust- 
free. The radial profiles in different bins of M* (right panel) show 

7 The difference between A RLR and ( A v ) ga,axy is -0.03 ± 0.06, while 
between A RLR and it i s _Q.0 ±0.1. 
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Fig. 15. Radial profiles of A v as a function of Hubble type (left), and in seven bins of galaxy stellar mass (right). These bins are log M*(M 0 ) = 
9.1-9.6, 9.6-10.1, 10.1-10.6, 10.6-10.9, 10.9-11.2, 11.2-11.5, 11.5—11.8. Symbols and colors are as Fig. 6. These results are obtained with 
base GMe. 


a similar behavior to that with morphology. Except for the most 
massive bins, shifted to lower extinction values, all other mass- 
binned Ay(R) profiles are similar. 

5.4.3. Radial gradients 

Ay gradients are shown in Fig. 16, which is formatted as Figs. 7, 
10 and 13. As for the previous properties, results for bases GMe 
and CBe are nearly indistinguishable, as illustrated by the over¬ 
lap of circles and stars in the left panel. V,„Ay and V our Ay show 
similar behavior with morphology, although the inner gradient is 
always higher than the outer one. In Ellipticals the gradient of 
Ay exists only in the central region. With the exception of Sd 
galaxies, spirals have V,„ Ay -0.25 mag/HLR. 

On average, V,„Ay gets stronger with increasing M* up 
to 10 n M o (Fig. 16, right) and weakens towards higher mass, 
spheroid dominated systems. The dispersion with respect to the 
mass-binned relation (traced by the black crosses) is large, and 
not clearly related to morphology (coded by the colored circles). 

As a whole, and despite the general trends summarized 
above, of the four properties analyzed in this section. Ay is the 
one for which tendencies with Hubble type and stellar mass 
are less clear. A major reason for this is that, unlike for /r*, 
(log age)L, and (logZ*)^, Ay estimates are sensitive to inclina¬ 
tion effects. This is explained next. 

5.5. Effect of inclination on the radial profiles 

An implicit hypothesis throughout the analysis presented so far 
is that galaxy inclination does not affect our estimates of the stel¬ 
lar population properties and their radial distributions. One ex¬ 
pects this assumption to break down in the case of Ay , which 
should increase from face on to edge on galaxies, although it is 
not unreasonable to conceive that inclination effects propagate 
to the spectral synthesis-based estimates of stellar mass surface 
densities, mean ages, and metallicities. It is therefore relevant to 
evaluate if and how inclination affects our results. 

In order to do so, we have divided the 300 galaxies in three 
subsamples on the basis of the b/a ratio (minor to major isopho- 
tal axes), as measured in SDSS /(-band images. The three sub¬ 
samples, each containing 100 galaxies, cover (i) b/a < 0.39, 


edge on, (ii) 0.39 < b/a < 0.63, and (Hi) b/a > 0.63, face on. 
Galaxies in each sub-sample were grouped by Hubble type, and 
their radial profiles of log/;*, (log age)L, (log Z*) jM , and Ay av¬ 
eraged as previously done for the whole sample in the left panels 
of Figs. 6, 9, 12, and 15. 

Fig. 17 shows the resulting stacked profiles of log /./», 
(log age)i, (log Z*) M , and A v . Solid, dashed and dotted lines 
show profiles for the “face-on" (b/a > 0.63), intermediate in¬ 
clination (0.39 < b/a < 0.63), and “edge-on” ( b/a < 0.39) sam¬ 
ples respectively, and each column is for one of the seven Hub¬ 
ble types used throughout the paper. Average profiles were only 
computed for morphology-inclination bins containing at least 4 
galaxies. 

Stellar mass surface density, age, and metallicity profiles 
show a negligible variation among the b/a -based subsamples. 
This result indicates that inclination does not affect the estimates 
of these properties in any systematic way. Any difference is at a 
level not significant insofar as this paper is concerned. The ex¬ 
ception is the (log Z*) M profiles for “edge-on” Sc’s, which differ 
substantially from the profiles of less inclined Sc’s. It so hap¬ 
pens, however, that the sub-group of b/a < 0.39 Sc’s has a mean 
stellar mass 0.4 dex lower than other Sc’s, which could explain 
their lower metallicities without implying inclination effects. 

The one property which does vary systematically with b/a 
is A v , and it does so in the expected sense: Spirals with lower 
b/a have larger extinction. This is particularly evident in Sb’s. 
This dependence hinders the interpretation of the stacking results 
presented in §5.4, and explains why no clean tendencies of the 
Ay values and profiles with morphology and stellar mass were 
identified. 

6. Discussion 

This section is divided into four main parts. First we summa¬ 
rize our results in the context of related studies. We then discuss 
our findings in the context of the growth of galaxies - theoreti¬ 
cal expectations, high redshift observations, and previous results 
of inside-out growth for CALIFA galaxies. In the third part we 
explore what the results tell us about the quenching of star for¬ 
mation in galaxies. Finally, we discuss the theoretical predictions 
for the radial variations of age and metallicity in early types and 
in spirals from different simulations of galaxy formation. We 
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Fig. 16. Left: As Fig. 7 but for Ay. Right: The inner gradient as a function of galaxy stellar mass. Symbols and colors are as in Fig. 7. 



radial distance (HLR) 


Fig. 17. Radial profiles of logyu*, (log age) L , (log Z*) M , and A v (from the upper to the bottom panels) for the different Hubble types (from E to Sd 
from left to right panel) and three different bins of the ratio of minor to major photometric radius of the galaxy: solid line (0.63 < b/a, face on), 
dashed line (0.39 < b/a < 0.63), and dotted line (b/a < 0.39, edge on). 

compare our results for variations of the radial structure in the 6.1. Age and metallicity of bulges and disks 
inner (R< 1 HLR) and outer (1 <R< 3 HLR) parts with other 

observational results in the literature. The analysis of SDSS data has generated a general knowledge of 

how the age and metallicity of galaxies change with M*, color, 
or concentration index (e.g Kauffmann et al. 2003; Gallazzi et al. 
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2005; Mateus et al. 2006). These studies have confirmed that, in 
general, early type galaxies are old and metal rich, while late 
type galaxies are younger and more metal poor. Numerous (sin¬ 
gle spectra or longslit) studies have reported also that ellipticals 
are metal rich, and have a range of ages, 2-10 Gyr, that depend 
on stellar velocity dispersion (e.g. Trager et al. 2000; Thomas 
et al. 2005; Sanchez-Blazquez et al. 2006; Graves et al. 2009; 
Johansson et al. 2012). 

Our spatially resolved analysis introduces a significant im¬ 
provement in the study of the structure of galaxies. For exam¬ 
ple, we compute ages and metallicities of bulges in disk galaxies 
and compare them with elliptical galaxies in a systematic way, 
avoiding problems derived from the lack of spatial resolution 
that some of the previous studies have. 

We compute the luminosity-weighted and the mass-weighted 
age and metallicity: (i) in the central part of galaxies (values av¬ 
eraged within 0.25 HLR) as representative of the stellar popula¬ 
tion properties of bulges and central core of ellipticals; and (ii) 
in the outer part of galaxies, values averaged in a ring at 1.5+0.1 
HLR, as representative of disks. Fig. 18 plots the individual re¬ 
sults as small dots; large dots represent the average values for 
each (color-coded) Hubble type, and the error bars show the dis¬ 
persion. While (log age)i gives information about the globally 
‘averaged’ star formation history, (log age )m informs when most 
of the stellar mass was formed. 

Fig. 18 shows that the bulges of Sa-Sb and the cores of E- 
S0 formed at a similar epoch; they are very old (~ 10 Gyr) and 
metal rich (>1 Z 0 ). Thus, they probably formed by common pro¬ 
cesses, that occurred rapidly and early on. However, the bulges 
in Sc-Sd galaxies (shown as the two darkest shade of blue) are 
younger and have a more extended star formation history (both 
(log age )m and (log age) l are smaller), and have lower stellar 
metallicities. Thus, Sc-Sd galaxies formed in later times and/or 
by different processes. 

Many bulges in late type spirals are in fact pseudo-bulges. 
Unlike true bulges, pseudo-bulges are thought to grow by secular 
processes, where material from the disk (by the loss of angular 
momentum) is the main source of star formation in the central 
parts of these galaxies (e.g. Kormendy & Kennicutt 2004). We 
may see this effect at work in Fig. 9 and Fig. 12, as a flattening of 
the radial profiles of (log age)i and (logZ*) M , and the positive 
(log age)L gradient in the core of Sc galaxies. Some effects of 
the secular processes due to the disk may also be present in the 
bulges of Sa-Sb. For example. Fig. 18 shows that bulges of Sa- 
Sb have (logage)i~6 Gyr, younger than the 10 Gyr epoch of 
formation derived from (log age)M~, and this may be understood 
if some disk stars are rearranged into the bulges or if dissipation 
processes bring gas triggering new star formation in the center. 

Fig. 18 also shows that disks are younger and more metal 
poor than bulges. Both (logage)i and (log age )m are lower in 
disks than in their respective bulges, indicating that disks formed 
later than bulges, and that star formation continues for a longer 
time in disks that in bulges, probably as a consequence of a con¬ 
tinuing availability of gas in disks (Roberts & Haynes 1994). 
This indicates a general scenario of inside-out formation. 

6.2. Inside-out growth of spheroids and spirals 

6.2.1. Theoretical expectations and recent results from high 
redshift galaxies 

Models of galaxy formation predict a common inside-out view 
for the mass assembly in galaxies (e.g. Kauffmann et al. 1993; 
Aumer & White 2013). First, the bulge formed at high redshift; 


then, the disk was built around the bulge in the case of spirals. 
In the case of ellipticals, the central core formed at z>2, and the 
envelope grew later through minor mergers (Oser et al. 2010; 
Hilz et al. 2013, e.g). Observational evidences come from the 
significant size evolution in early type galaxies (ETG), that grow 
as size oc M* (van Dokkum et al. 2010; Patel et al. 2013). 

More recently, van Dokkum et al. (2014) find evidence 
against the inside-out formation scenario for spirals. For a sam¬ 
ple of MW-like spirals at redshift z=2.5, they estimate the depen¬ 
dence of the radius with M*, and find that their size-M* relation 
is similar to the size-M* of similar galaxies at z=0. They con¬ 
clude that the mass growth took place in a fairly uniform way, 
with the galaxies increasing their mass at all radii, thus, their 
R e ff barely grows. These results seem to be supported by numer¬ 
ical simulations by Elmegreen et al. (2008), that find that bulges 
can be formed by migration of unstable disks. Other observa¬ 
tional evidence come from the detection of clumpy star forming 
disks in galaxies at z~2 (Genzel et al. 2008; Forster Schreiber 
et al. 2011), that may indicate an early build up of bulges by sec¬ 
ular evolution. Thus, studies at high redshift are providing new 
results that draw a complex landscape of galaxy build up. For ex¬ 
ample, Wuyts et al. (2011) also find clumpy disk star formation, 
but at the same time conclude that there is a Hubble sequence in 
place at least since z~2.5. On the other hand, there is other ev¬ 
idence that galaxies rapidly assemble inside-out at z= 1 (Nelson 
et al. 2012; Szomoru et al. 2010, 2012); while Hammer et al. 
(2005) find evidence that MW-like galaxies have rebuilt their 
disk at z < 1 in a major merger epoch that drastically reshapes 
their bulges and disks, and is consistent with earlier cumplier 
evolution. 

In summary, there is mounting evidence of the major pro¬ 
cesses responsible for the assembly and shaping of galaxies at 
different epochs, and these are complemented with a variety of 
processes that modify the inside-out formation scenario: stellar 
migration, bar induced gas inflows, gas-rich minor merger, angu¬ 
lar momentum loss due to reorientation of the disk, infall of gas 
with misaligned angular momentum, etc (Aumer et al. 2014). 

6.2.2. CALIFA view of the inside-out growth of galaxies 

The results from our studies favor an inside-out growth of spi¬ 
rals. Perez et al. (2013) studied the stellar mass growth as a func¬ 
tion of the radius and cosmic time in galaxies with 10 10 < M* < 
5 X 10" M 0 , and showed that the nuclei grow faster than the 
inner 0.5 HLR, that, in turn, grow faster than the outer 1 HLR. 
This conclusion is supported by the stellar age radial profiles pre¬ 
sented in Gonzalez Delgado et al. (2014c), and confirmed here in 
Fig. 10 for most spirals and spheroidals. Further support comes 
from the ratio HMR/HLR (Fig. 3), a good probe of the spatial 
variation of the star formation history in galaxies (Gonzalez Del¬ 
gado et al. 2014c). This ratio is lower than 1 (Fig. 3), a fingerprint 
of the inside-out growth found by Perez et al. (2013). 

Fig. 19 shows how the radial profiles of (log age)M decrease 
outwards for all the Hubble types. Most of the stellar mass in 
the center has been formed 10 Gyr ago or earlier (z>2). But at 
1.5 HLR, (log age)M ranges from 7 Gyr (z~ I ) in E-S0 to 4.5 
Gyr (z~0.4) in Sbc, suggesting that, both early type and MW- 
like, galaxies have continued accreting or forming in-situ stars 
in their disks until more recent times, thus supporting the inside- 
out scenario in these galaxies. 

This trend, however, changes beyond 1.5-2 HLR, where 
(log age)M and (log age )l flatten. This may be interpreted as 
indicating that the mass was formed in a more uniformly dis¬ 
tributed manner across the outer disk, or that stellar migration 
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Fig. 18. Left: Mass-weighted age at the galaxy center ({log age) M bulge ) and at 1.5 HLR ((log age) M d,sk ) for the different Hubble types. Small dots 
are the individual radial points, while big colored dots represent mean values for each Hubble type, with the error bars indicating the dispersion in 
the mean. Middle: As in the left panel but for the light-weighted age ((log age) L ). Right: As in the left panel for (logZ*)^. The top horizontal axes 
in the left and middle panels show the redshift scale. The diagonal line in the three panels is the bisector. 



Fig. 19. Radial profiles (in units of HLR) of the mass weighted age, 
(log age)M bulge , obtained with GMe base. The results are stacked by 
morphological type as in Fig.7. 

shuffles inner born stars to the outer disk, washing out the inside- 
out formation signs so clearly seen in the inner 1.5 HLR. In the 
case of E-S0 this may be understood if beyond 2 HLR most of 
the stellar mass in the galaxies was already accreted at z < 1. 

6.3. Quenching 

Several mechanisms have been proposed to explain the shut¬ 
down of star formation in galaxies. Halo mass quenching is one 
of the most popular ones that explains the bimodal distribution 
of the properties of galaxies, and it is required to explain the 
green valley as a pathway towards quenching of star formation 
in early and late type galaxies (e.g. Schawinski et al. 2014). In 
this scheme, galaxies with a halo mass below a critical value 
(a few x 10 12 M g ) accrete cold gas conducive to star forma¬ 
tion. Above this critical mass, the infalling gas reaches the sound 
speed and cannot form new stars (e.g. Cattaneo et al. 2006; Dekel 
& Birnboim 2006). The dependence with environment and clus¬ 
tering strongly supports this quenching mechanism (e.g. Wein- 
mann et al. 2006; Peng et al. 2010). 

The differential dependence of the stellar mass surface den¬ 
sity (Fig. 5 and Fig. 6) with the galaxy stellar mass (a proxy of 
the halo mass) provides further evidence of the halo quenching 
(e.g. Behroozi et al. 2010). Estimating of the properties of the 


stellar populations in SDSS galaxies, Kauffmann et al. (2003) 
found that there is a critical mass (Af* = 3 x 10 10 M Q , equiva¬ 
lent to ~ 6 x 10 10 M g for our Salpeter IMF) below which log//, 
scales with M*. and above which log//* is independent of the 
galaxy stellar mass. Right panels of Fig. 5 and Fig. 6 support this 
scenario because the radial profiles of log //* scale with log Af*, 
and furthermore they do so all along their extent. Our results also 
show that log//* saturates at high Af*; because the high mass end 
of the distribution is dominated by early type galaxies (Sa-SO-E), 
this suggests that the spheroidal component plays a significant 
role in the quenching of star formation in high mass galaxies. 

The importance of morphology in the quenching of galaxies 
has also been reported in the literature (e.g. Bell 2008; Bell et al. 
2012; Barro et al. 2013; Pan et al. 2014; Woo et al. 2015). Martig 
et al. (2009) found that the dependence of quenching with mor¬ 
phology is a consequence of the bulge-building mechanism. The 
steep potential well induced by the formation of a large spheroid 
component results in the stabilization of the disk, that cuts the 
supply of the gas, preventing its fragmentation into bound, star 
forming clumps. Our results support this scenario, as it is ex¬ 
plained below, because the dependence of the SFH of galaxies 
with the morphology. 

If the halo mass is the main property responsible for quench¬ 
ing, we should expect that the radial structure of (log age)^ 
(both, the age values and the gradients) to change more with Af* 
than with Hubble type. On the contrary, if quenching is driven 
by morphology, galaxies of similar stellar mass would have very 
different (log age )/, structure depending on Hubble type. We ex¬ 
plore the relevance of morphology versus Af* in Fig. 20: age ra¬ 
dial profiles are shown as a function of Af* and of morphology, 
in four mass bins (log M* = l 1.5-11.2, 11.2-10.9, 10.9-10.6, 
10.6-10.1). Clearly, morphology is the main driver: it can ac¬ 
count for up to 0.75 dex change in age at a given mass (top pan¬ 
els); conversely, at a fixed morphology, mass accounts for less 
than 0.25 dex (bottom panels). Further, morphology accounts not 
only for changes in absolute values, but also for changes in the 
gradients at a given galaxy mass. 

This confirms the similar result obtained above with log//*, 
and it implies that galaxies of similar Af* (equivalent to have 
similar Af/ M / 0 ) and with a large spheroid have shutdown their 
star formation (outside their central core) earlier than galaxies 
of later morphology. These results indicate that the SFH and 
their radial variations are modulated primarily by galaxy mor¬ 
phology, and only secondarily by the galaxy mass, suggesting 
that the bulge formation has a relevant role in quenching the star 
formation in galaxies. 
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Fig. 20. Radial profiles of (log age)/, (upper panel) in four galaxy stellar mass bins. From left to right: log M*(M G ) = 11.2-11.5 (continuum line), 
10.9-11.2 (dashed line), 10.6-10.9 (dashed-point line), 10.1-10.6 (dotted line). In each panel, the average profile for each Hubble type is plotted 
if more than four galaxies have galaxy stellar mass in the log M* bin. Bottom: each panel shows the radial profile of each Hubble type averaged in 
each of the four log M*(M G ) bins. 


6.4. Radial structure of the stellar population properties in 
ETG and their relation with galaxy formation models 

6.4.1. Theoretical predictions from cosmological simulations 

Classical chemical evolution models of the formation of early 
type galaxies (ETG) are based in two possible scenarios: 1) dissi¬ 
pative formation, the well known monolithic collapse; and 2) the 
non-dissipative collapse. These scenarios produce very different 
radial gradients of ages and abundances, being very steep in the 
first case, with values of V[Fe/H]-0.5 to -1.0 [dex/dex] (Lar¬ 

son 1974, 1975; Carlberg 1984) 8 , but (almost) flat when there 
are pure stellar mergers. This second case may even erase a pre¬ 
viously existing radial gradient. 

The most recent cosmological simulations propose a two 
phase formation scenario for ETG’s in which the central core 
formed at z > 2, and the envelope grows after this through mi¬ 
nor mergers (e.g. Naab et al. 2009; Oser et al. 2012; Hilz et al. 
2012; Navarro-Gonzalez et al. 2013). Thus: 1) Galaxies assem¬ 
ble their mass through dissipative processes and star formation 
occurs in-situ. Starbursts formed at the center as a consequence, 
for example, of large major mergers or monolithic collapse. The 


8 The metallicity gradient measured in spheroids is traditionally calcu¬ 
lated as A[Fe/H]/A log r and expressed in [dex/dex] units. 


star formation is induced by cold flow of accretion or by gas-rich 
mergers. 2) Galaxies grow in size by mass assembly through the 
external accretion of satellites; ex-situ star formation formed by 
dry mergers of galaxies towards the central most massive one. 

Observationally, there is evidence of a significant size evolu¬ 
tion in ETGs. The growth of the galaxy size with supports 
this proposal. A transition region is expected between the in-situ 
central core of ETG and ex-situ outer regions. Since the central 
core of these ETG is enriched very quickly due to major merg¬ 
ers at high redshift (z > 2), and the satellites that are accreted 
are less metal rich than the central massive one, a negative ra¬ 
dial gradient of the metallicity is expected, even with a change 
of the slope in the transition region. Thus, values as V[Fe/H] = 
-0.5 [dex/dex] (Pipino et al. 2010) or V[Fe/H] = -0.3 [dex/dex] 
(Kawata & Gibson 2003) are predicted. 

However, the merger history may change an existing radial 
gradient: while dry major mergers can flatten the pre-existing 
gradient (Kobayashi 2004; Di Matteo et al. 2009; Hopkins et al. 
2009), dry minor mergers can steepen the metallicity gradient. 
Thus, Kobayashi (2004) SPH chemodynamical simulations of 
elliptical galaxies that include radiative cooling, star formation 
and feedback from SNII-Ia, and chemical enrichment, (but do 
not include kinematic feedback), found that the steep negative 
radial metallicity gradient, established during the initial starburst 
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at z > 2, can flatten significantly by later major-mergers in the 
galaxy. Following these simulations, the average gradient at the 
present time is V[Fe/H]= -0.3 [dex/dex], but it may decrease to 
a value of -0.2 [dex/dex] when major mergers appear. 

Beside the merger history, feedback can change the inner 
and outer metallicity gradients. Thus, a strong AGN feedback 
can stop the star formation in the central core of massive galax¬ 
ies, flattening the inner gradients. Feedback from in-situ star for¬ 
mation can alter the outer metallicity gradient. Also, the exis¬ 
tence of galactic winds may modify the composition of the ISM 
in a galaxy. Hirschmann et al. (2015) performed cosmological 
simulations that include an empirical model for the momentum 
driven galactic winds, to investigate the dependence of the age 
and metallicity outer gradients with metal cooling and galactic 
winds, (in principle required to explain the mass-metallicity rela¬ 
tion, MZR). These simulations including winds predict V[Fe/H] 
= -0.33 [dex/dex], steeper than the simulations without winds 
that predict V[Fe/H] = -0.11 [dex/dex]. The main explanation 
is that in wind models the stars accreted are of lower metallicity 
than in the simulations with no winds. In both cases, however, 
they predict a positive age gradient of ~ 0.03-0.04 [dex/dex]. 

6.4.2. Implications from this work and comparison with other 
results from the literature 

Following our own results in this work, E and SO have formed 
their central core at similar cosmic time since they have sim¬ 
ilar central ages (see Fig. 18). Further they must have formed 
through similar processes since their radial profiles of log/r*, 
(log age)L, and (log Z*) ;W are remarkably similar. They both 
show small but negative (logage) l, and (logZ*)^ gradients in 
the central core. In the central 1 HLR, E and SO in our sam¬ 
ple have V(logZ*)M~ -0.1 [dex/dex] (std = 0.15) 9 . (Sligthly 
steeper, V(log Z*)m= -0.2 [dex/dex], when CBe models are 
used.) These are within the range of values found in other stud¬ 
ies based on long-slit or IFS data up to one effective radius (e.g. 
Mehlert et al. 2003; Sanchez-Blazquez et al. 2007; Annibali et al. 
2007; Rawle et al. 2008; Spolaor et al. 2010; Kuntschner et al. 
2010). However they are shallow compared with theoretical ex¬ 
pectations if minor mergers are relevant in growing up the central 
core of E and SO galaxies. This may indicate that major mergers 
are more likely the main process building the central regions (up 
1 HLR) of ETGs. 

Between 1 and 3 HLR, the radial profile of (logZ*)^ is of 
similar slope or slightly shallower than in the inner 1 HLR. We 
do not find any evidence of a transition region where the metal¬ 
licity radial profile steepens to metallicities below solar. If the 
1 to 3 HLR envelope of ETG had grown through the accretion 
of low mass satellites a steepening of metallicity would be ex¬ 
pected, as explained before, because the mass-metallicity rela¬ 
tion implies that low mass satellites would be of low metallicity. 
In our results there is no evidence either of an inversion of the 
age radial profile toward older ages beyond 1-2 HLR, as ex¬ 
pected if these satellites were formed very early on like the core 
of E and SO (see Fig. 9 and Fig. 18). These results are in contrast 
with recent ones by Greene et al. (2012, 2013): for a sample of 
~30 early type galaxies they find at 2 R,,// an old (10 Gyr) stel¬ 
lar population with [Fe/H]-0.5, and interpret this as the stellar 

outskirts of these galaxies being built up through minor merg¬ 
ers. Also Coccato et al. (2010); La Barbera et al. (2012); Montes 


9 Our gradients, that are measured in a linear scale, are converted here 
to a logarithmic scale to be compared with predictions front simulations 
and other works in the literature. 


et al. (2014) observing a few massive ellipticals have reported a 
decline of the metallicity to under solar in an old stellar popula¬ 
tion in their outskirts (> 10 Re//) suggesting that these galaxies 
are formed in two phases, the central core through major merg¬ 
ers, and through minor mergers farther out. However, other re¬ 
cent works show examples of ETGs with an old and metal rich 
stellar population and a very shallow metallicity gradient up to 3 
R e // (Trujillo et al. 2014), in contrast with the results by Greene 
et al. (2012, 2013). 

Our results do not support the minor merger scenario for the 
size growth of ETGs. Thus, the ages, (log age)^~ 9.7 (yr), and 
metallicity, (log Z*)m~ Z 0 , at 1 - 3 HLR, and the shallow metal¬ 
licity gradient, V/logZ*)^-0.1 [dex/dex], that we obtain are 

more consistent with the growth of the 1-3 HLR envelope of 
ETGs through major mergers. 

Other interesting result reported in the literature that can be 
compared with ours is the correlation between the metallicity 
gradient and the galaxy mass (or stellar velocity dispersion, cr*) 
found for E and SO. Spolaor et al. (2010) have found that the rela¬ 
tion between <x* and the metallicity gradient shows two regimes 
of behavior: (i) for galaxies with log cr* < 2.2 km s', the metal¬ 
licity gradient steepens with cr*; (ii) galaxies with log cr* > 2.2 
km s 1 (the most massive ellipticals), have a metallicity gradient 
that does not correlate with cr* (or galaxy mass), with a mean 
value ~ -0.15 [dex/dex]. On the other hand, Pastorello et al. 
(2014) derive the [Z/H\ gradient in the outer 1 - 2.5 R<,// of a 
sample of ellipticals, and they find that the gradient covers a wide 

range of values, from negative very steep (-2) to flat or even 

positive; these values correlate with the galaxy stellar mass and 
stellar velocity dispersion, with the galaxies of lower M* (or cr*) 
having the steeper metallicity gradient. However, the most mas¬ 
sive galaxies exhibit the flattest gradients and an average value 
of V(log Z*)m~ -0.2 [dex/dex] (std = 0.38) for galaxies with 
M* > 10 11 M 0 . Both works show a significant scatter in the 
V(log Z*)— M* relation for M* > 10 11 M 0 , and reasonable 
doubts of the existence of the correlation for high mass ellipti¬ 
cals. 

Our results (Fig. 13) indicate that there is no correlation be¬ 
tween the metallicity gradient and M* for the CALIFA early 
type galaxies (E and SO). Even so, our results are compatible 
with Spolaor et al. (2010) and Pastorello et al. (2014), because 
E and SO in our sample are all above 10 11 M 0 (and cr* >100 
km s' 1 ), for which no correlation is found between the metallic¬ 
ity gradient and M* in Spolaor et al. (2010) or Pastorello et al. 
(2014). We find that CALIFA ellipticals have a shallow gradi¬ 
ent. This behavior is also in agreement with Hirschmann et al. 
(2013) simulations and their interpretation of the lack of corre¬ 
lation between the metallicity gradient and M* in massive ellip¬ 
ticals. Following these authors, massive galaxies accrete higher 
mass satellites, and because of their deeper potential well they 
retain their own gas against stellar winds, producing a shallower 
metallicity gradient in the outer regions of massive ellipticals. 

Because the metallicity at 2-3 HLR in E and SO are similar 
to the metallicity in the bulge of early spirals, and the stars at 
these distances are as old as the bulges of Sa-Sb galaxies (see 
Fig. 9), the 1-3 HLR envelope of early type galaxies might have 
built from the centers of early type spirals. In summary, the nega¬ 
tive but shallow gradients of the metallicity and ages suggest that 
massive (M* > 10 11 M 0 ) early type galaxies built their inner 3 
HLR through mergers with massive and metal rich spirals. 
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6.5. Radial structure of the stellar population properties in 
spirals and their relation with galaxy formation models 

New insights on the structure of the Milky Way disk, in partic¬ 
ular through the measurements of chemical abundances of large 
sample of stars, are provided by the spectroscopic surveys under¬ 
taken in recent times (e.g., SEGUE, RAVE, Gaia-ESO survey, 
HERMES, APOGEE, LAMOST, etc; Yanny et al. 2009; Stein- 
metz et al. 2006; Gilmore et al. 2012; Zucker et al. 2012; Majew- 
ski et al. 2010; Zhao et al. 2012). RAVE (Radial Velocity Exper¬ 
iment) (Steinmetz et al. 2006; Boeche et al. 2014) is studying the 
radial and vertical chemical gradients using a very large sample 
of dwarf stars. Close to the Galactic plane, RAVE shows a neg¬ 
ative radial gradient of Fe abundance, -0.054 dex/kpc 10 , that 
becomes flatter or even positive when measured above the disk. 
So, the [Fe/H] gradient ranges from -0.039 to +0.047 dex/kpc 
when measured at heights 0.4 - 0.8 kpc, or 1.2 - 2 kpc above the 
Galactic plane, respectively. 

The radial gradient of abundances in the different regions of a 
spiral galaxy are important because they are directly related with 
the formation process. Obviously, not all scenarios of disk/spiral 
formation are valid, since it is necessary that they produce a ra¬ 
dial gradient of abundances in the disk but not in the halo, as 
observed in the MW and in M31. Thus, the formation of the 
halo from different fragments or minor mergers with very short 
free fall times does not create a radial gradient but a dispersion 
of abundances, and therefore it was early concluded that the MW 
halo may be formed from mergers or from the accretion of low 
mass galaxies (or part of them). However, disks are more likely 
formed from a single cloud falling on and from inside-out. 

6.5.1. Theoretical predictions from "classical" chemical 
evolution models 

Most classical chemical evolution models claim that infall of gas 
with a radial dependence, implying an inside-out scenario for the 
disk formation, is essential to reproduce the observed radial gra¬ 
dient of abundances. The key ingredient is the dependence of the 
disk infall time scale with the radial distance, that makes the gas 
to accumulate faster in the inner disk. Since the SFR depends 
on the gas density, these assumptions produce a radial depen¬ 
dence of the star formation rate and a negative radial metallic- 
ity gradient (Ferrini et al. 1994; Molla et al. 1996; Chiappini 
et al. 2001; Molla & Diaz 2005). Thus Molla et al. (1996) give a 
value -0.08 for a MW like galaxy, reproducing the value found 
by Shaver et al. (1983) and other Hu regions studies. Chiappini 
et al. (2001) models predict a gradient of -0.04 dex/kpc for a 
MW-like galaxy. In fact, as theoretical equations show (Goetz & 
Koeppen 1992), the radial gradient of abundances appears in the 
disk when there is an adequate ratio between star formation rate 
to infall rate. It also implies, therefore, that a dependence of the 
radial gradient on the morphological type of galaxies may ex¬ 
ist. Molla et al. (1996) models already predicted radial gradients 
for galaxies of different morphological types, with values in the 
range -0.025 (for a M31-like galaxy) to -0.183 dex kpc 1 (for 
a late type galaxy like NGC 300). More recent works (Molla & 
Diaz 2005) calculate models where the infall rate was a function 
of the mass distribution (or rotation curve) of the galaxy, assum¬ 
ing a stronger radial dependence of the infall timescales than in 
Chiappini et al. (2001). Moreover Molla & Diaz (2005) models 
also depend on an efficiency factor to condense the molecular 


111 The metallicity gradient in disks is traditionally calculated as 
A[Fe/H]/Ar, and expressed in dex/kpc. 


gas, and to convert the gas reservoir into stars. The metallic¬ 
ity gradients range -0.02 to -0.15 dex kpc 1 with flat gradients 
for galaxies with the largest efficiency factor, or the most mas¬ 
sive ones, although in the extreme end the low mass and low¬ 
est efficiencies models also show flat radial distributions. Thus, 
the steepest gradients appear in the intermediate mass or inter¬ 
mediate type galaxies. However, there is no dependence on the 
morphological type when the gradient is normalized to a char¬ 
acteristic value, such as the effective radius as recent results by 
CALIFA have found based on HII regions abundances (Sanchez 
etal. 2013). 

6.5.2. Theoretical predictions from cosmological simulations 

Recently, hydrodynamical cosmological simulations have pro¬ 
vided evidences in support of the imposed inside-out disk growth 
scenario adopted within the "classical" chemical evolution mod¬ 
els. Like spheroidals, spirals are formed in two phases. In the 
first phase the bulge formed in a similar way as the core of E- 
S0. In the second phase, the disk grows by star formation in-situ 
from the infalling gas (Kauffmann et al. 1993; Aumer & White 
2013). Metal poor gas with higher angular momentum at lower 
redshifts is turned into stars at larger radii. Negative radial metal¬ 
licity gradients are expected, as the classical models predict. This 
assumption is a natural outcome of the mass, momentum, and 
energy conservation laws, imposed in the simulations of disks 
in a cosmological context (Brook et al. 2011, 2012; Few et al. 
2012; Pilkington et al. 2012a,b; Gibson et al. 2013). 

Pilkington et al. (2012a) have examined a set of 25 simu¬ 
lations, from several groups, using different codes and initial 
conditions (Stinson et al. 2010; Rahimi et al. 2011; Kobayashi 
& Nakasato 2011; Few et al. 2012) to predict the present-day 
metallicity gradient in MW-like galaxies and its evolution. Al¬ 
though the evolution of the simulated metallicity gradients de¬ 
pends strongly on the choice of the sub-grid physics employed, 
most of the simulated galaxies tend to a similar present-day gra¬ 
dient of-0.05dex kpc ', in agreement with the Chiappini et al. 

(2001) and Molla & Diaz (2005) models for normal galaxies as 
the MW. 

6.5.3. Implications from this work 

Our findings show that spiral galaxies (excluding Sd) have nega¬ 
tive radial gradients as indicative of the inside-out growth of the 
disk (see Fig. 12). The average V ouf (log Z*)m for spirals (exclud¬ 
ing Sd and later type) is -0.08 dex/HLR or -0.02 dex/kpc. These 
values are compatible with the results obtained by Sanchez- 
Blazquez et al. (2014), that have already derived the metallicity 
gradients for 62 CALIFA face-on spiral galaxies to study the ef¬ 
fect of bars on the properties of the stellar populations. For these 
galaxies, they find a metallicity gradient of -0.025 dex/kpc (std 
= 0.05), equal (-0.027 dex/kpc) to the gradient that we derive 
for the same group of galaxies. 

In order to compare with our results, the RaDES simulated 
galaxies (Few et al. 2012), 19 galaxies of the Pilkington et al. 
(2012a) sample, have been analyzed in a similar way as we have 
done here, i.e. measuring the gradients in a similar way and 
using the HLR values of the simulated galaxies (Ruiz-Lara in 
prep., and Ruiz-Lara et al. private communication). The simu¬ 
lated galaxies analyzed in this work cover a narrow range of 
morphologies, mainly Sbc-Sc. Therefore, these results cannot 
be extrapolated to the full work presented here, but they are 
representative of the state-of-art of cosmological simulations of 
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disk galaxies, and can be used to compare them to similar disks 
from our observations. Mock B-band images are used to derive 
the HLR and to perform a bulge-disk decomposition used as a 
proxy for the morphology. Metallicities are calculated for disk 
particles using Eq.(2) and the gradient is derived between 0-1 
HLR and also between 1 and 2 HLR. The stellar metallicity and 
age gradients of the simulated galaxies are compatible with the 
results presented here. Keeping in mind that the morphological 
range covered by these simulations is rather narrow and that they 
use B/D as a proxy for a morphological classification, the re¬ 
sults show a slight dependence of V ou , (log Z+)m with B/D ratio, 
with a steeper slope for B/D= 1 (Sbc galaxies) for which V OM , 
(logZ*) M -0.1 dex/HLR. Later type spirals have a flatter gra¬ 
dient of V OM , (log Z*),v;-0.046 dex/HLR. These results go in 

line with those found here, namely, that the metallicity radial 
gradient of spirals shows a dependence on morphology, with the 
steepest gradient found in the intermediate Sb-Sbc spirals (see 
Lig. 12 and Lig. 13). However, a larger set of cosmological sim¬ 
ulations is required covering from early Sa to late Sd, and a large 
range of galaxy masses (from 10 9 to 10 11 M G ), in order to con¬ 
firm the general trend found here. On the other hand, these re¬ 
sults indicate that the feedback recipes used in these simulations 
are able to recover realistic galaxies with small bulges and are 
fully in agreement with the work presented here. 

Lurthermore, our results are also compatible with classical 
chemical models, and certainly, the CALILA Sb-Sbc galaxies 
have stellar metallicity gradients (-0.025 dex/kpc) in the range 
observed in the MW disk, but somewhat shallower than the 
[Le/H] gradient measured in the Galactic disk. However, it is 
necessary to take into account that the gradient usually given in 
the literature is obtained for young stars or HII regions, while 
here it is an average value obtained for all stellar populations ex¬ 
isting in the studied galaxy or region. Besides that, the number 
of objects is increased compared with the old studies and, more 
important, all of them have been self-consistently analyzed using 
the same reduction technique and spectral models. 

In any case our results favor an inside-out growth of spirals. 
This conclusion is supported by the stellar age radial profiles 
presented here: the age decreases outwards for all Hubble types 
studied 11 . Beyond 1.5-2 HLR the radial distribution of ages 
flattens, suggesting that the mass forms more uniformly in those 
regions, or that the stellar mixing brings stars born in the inner 
disk to the outskirts. This last possibility has been recently in¬ 
vestigated by Minchev et al. (2013, 2014), who have performed 
N-body hydrodynamical models with the chemical evolution im¬ 
plementation (Minchev et al. 2013, 2014). They have simulated 
MW-like galaxies with the aim to investigate whether the Galac¬ 
tic disk can be understood as a single structure with kinematic 
and chemical features that are continuously distributed, being 
the thin and thick disks two extreme cases of these structures. 
Lurthermore, they investigate the effect of stellar migration and 
kinematic heating in the scatter of the age-metallicity relation, 
and how it changes with the Galactic radius. In fact, an increase 
of the scatter in the age-metallicity relation and a flattening of 
the stellar metallicity gradient is produced by the stellar radial 
migration, that causes a radial mixing in the older stellar popula¬ 
tion, creating the appearance of a flatter gradient in early times, 
and leading to a decoupling of the stelar population from their 
birth interstellar medium (Roskar et al. 2008). These results also 
indicate that even though radial mixing has a significant effect in 
flattening the metallicity gradient, it can not destroy it. 


11 Sd galaxies, however, show a much flatter age gradient 


7. Summary and conclusions 

We have analyzed the stellar population properties of 300 galax¬ 
ies, observed by CALILA with the V500 and V1200 gratings 
and ILU PPak at the 3.5m telescope of Calar Alto, to investi¬ 
gate the trends in the stellar populations properties with radial 
distance as a function of Hubble type and galaxy stellar mass. 
The sample includes ellipticals, SO and spirals from early (Sa- 
Sb) to late types (Sc-Sd). They cover a stellar mass range from 
0.7xl0 9 to 7xlO u M 0 if Salpeter IMP is assumed, and a factor 
1.78 (0.25 dex) lower for a Chabrier IMP. A full spectral fitting 
analysis was performed using the starlight code and a com¬ 
bination of SSP spectra from Gonzalez Delgado et al. (2005), 
Vazdekis et al. (2010), or Chariot & Bruzual (2007, private com¬ 
munication). Our pipeline pycAsso is used to process the spectral 
fitting results to produce present day maps of the spatial dis¬ 
tribution of the stellar population properties. Lor each galaxy, 
these maps are azimuthally averaged to produce radial profiles 
(in units of the half light radius, HLR: c/( 0 ) of the stellar mass 
surface density (log/r*), stellar ages (light weighted, (log age)i, 
and mass weighted, (log age ) M ), metallicity ((logZ*) M ), and ex¬ 
tinction (Ay). The radial profiles are stacked as a function of 
Hubble type and of galaxy mass. Radial gradients of these prop¬ 
erties measured within the inner 1 HLR and between 1 and 2 
HLR are also obtained. 

Our main results are: 

1. Spatially averaged vs. integrated galaxy properties: the 
metallicity, (log Z*) M , galaxy-wide spatially averaged 
matches the metallicity obtained from the integrated spec¬ 
trum, and the metallicity at R=1 HLR. This result is equiv¬ 
alent to that obtained for the other stellar population prop¬ 
erties, log/.(*, (log age)i, and Ay, as reported by Gonzalez 
Delgado et al. (2014c,b), proving that effective radii are in¬ 
deed effective. 

2. Mass weighted size: We confirm our earlier finding 
(Gonzalez Delgado et al. 2014c) that galaxies are more com¬ 
pact in mass than in light by -20%. The HMR/HLR ratio 
shows a dual distribution with Hubble type, that breaks in the 
Sb-Sbc, the galaxies with the smaller HMR/HLR. This ratio 
also shows a dual dependence with M*: it decreases with 
increasing mass for disk galaxies, and becomes almost con¬ 
stant in spheroidal galaxies. These results are a signpost of 
the inside-out growth previously found by Perez et al. (2013). 

3. Stellar mass surface density: log/r*(R) shows declining pro¬ 

files that scale with morphology and with M*; this behavior 
is preserved at any given distance. At constant M*, log p*(R) 
is higher in early type than in late type spirals. E’s and SO’s 
show equal log p+Ut) profiles, independently of M*. The in¬ 
ner gradient, V,„ logcorrelates with Hubble type. The 
negative gradients steepen from late type spirals to spheroids, 
as well as with galaxy total mass in galaxies with M* < 10 11 
M q . At a constant M*, V,„ log/r* steepens with morphology, 
with E’s and SO’s having the steepest gradients. These results 
indicate that morphology, and not only M*, plays a relevant 
role in defining and the relation. 

4. Stellar ages: (logage)i(R) shows declining profiles that 
scale with morphology; this behavior is preserved at any 
given distance. Early type spirals are always older than late 
spirals. E’s and SO’s, although older than spirals, have both 
similar (log flge)i(R) profiles, indicating that these galax- 
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ies have similar star formation histories. The more massive 
galaxies are also the older ones; this “downsizing" behav¬ 
ior is always preserved at any given distance. The nega¬ 
tive V,„(log age)i depends on Hubble type in different ways: 
steeper from E and SO to Sbc, and shallower from Sbc to 
Sd. Thus, Milky Way like galaxies have the steepest age 
gradient. A V,„(log age)i,—M i , relation exists, increasing the 
gradient from the low mass galaxies (which have roughly 
flat profiles) up to about lO n M 0 , at this point the trend 
reverses and V,„(log age)L decreases with increasing M*. 
However, the dispersion in the V,„(log age')L-M i , relation 
and (log age)L HLR -M+ is significant and it is strongly re¬ 
lated with the morphology. Even more, the dispersion of the 
(log age)i(R) profiles of galaxies of equal mass is signifi¬ 
cant and larger than between the (logage)i(R) profiles of 
galaxies of different M* but the same Hubble type. Thus, the 
SFHs and their radial variations are modulated primarily by 
the Hubble type, with mass playing a secondary role. 

5. Stellar metallicity: (logZ*)M(R) shows mildly decreasing 
profiles for most Hubble types, except Sd’s that show little, if 
any, radial dependence. Milky Way like galaxies (Sbc) stand 
out as the ones with the steepest radial profiles, (log Z*)m(R) 
scales with M* in a similar way as it does with morphology. 
This can be understood as a consequence of the global mass 
metallicity relation -a primary dependence of the metallicity 
with M*. The metallicity gradients are negative but shallow 
on average, with V„,(log Z*),v/ ~ -0.1 dex/HLR, and show 
a small dependence with M* up to M* ~ 10 11 M 0 , steep¬ 
ening with increasing mass. Above 10 11 M 0 (Sa’s, SO’s and 
E’s) they have similar metallicity gradient. The dispersion 
in the V,„(log Z„) M - M* relation is significant and a trend 
with morphology is seen, in the sense that, for a given mass, 
intermediate type spirals are the ones with steeper gradients. 

6. Stellar extinction: All the galaxies show Ay(R) declining 

profiles, but do not have a clear trend with morphology or 
with galaxy mass. Most spirals show similar radial varia¬ 
tions, and similar average extinction. Ay ~ 0.2 mag at the 
disk and up to 0.6 mag in the center, with the inner gradi¬ 
ent V,-„A v-0.25 mag/HLR. However, Sd galaxies show 

a shallow central gradient. E and SO also show a negative 
gradient in the inner 1 HLR (shallower than in early type 
spirals), but out of the core they are almost dust free. On av¬ 
erage, v,-„Av gets stronger with increasing M* up to 10 11 M 0 , 
and weakens towards higher mass. However, the dispersion 
with respect to the binned mass relation is not related with 
Hubble type. A major reason for this is that Ay(R) profiles 
are sensitive to inclination effects, unlike (log age)i, or 
(log Z*) vf. Thus, spirals with larger inclination have larger 
extinction. This is particularly evident in Sb’s, that have the 
largest A y LR and the steepest central gradient. 

From these results, we conclude: 

- Evidence in favor of the inside-out growth of galaxies is 
found in the negative radial stellar age gradients. Metallicity 
gradients and the fact that galaxies are more compact in mass 
than in light also support this scenario. On the other hand, the 
flattening of the (log age)r(R) profiles beyond 1.5-2 HLR 
may be interpreted as indicative that the mass was formed in 
a more uniformly distributed manner across the outer disk of 
spirals. In the case of E’s and SO’s this may be understood if 


beyond 2 HLR most of the stellar mass was accreted at z< 
1 . 

- The mean stellar ages of disks and bulges are correlated, with 
disks covering a large range of ages, and late type spirals 
hosting the younger disks. The bulges of SO and early type 
spirals are old and metal rich as the core of E’s. They formed 
by similar processes, through mergers. Later type spirals, 
however, have younger bulges, and larger contribution from 
secular evolution are expected. Disks are younger and more 
metal poor than bulges, as an indicative of the inside-out for¬ 
mation scenario of these galaxies. 

- SO’s in this sample (all are massive galaxies), act as a transi¬ 
tion class between E’s and spirals, with //»(R ), (log age)i(R), 
and Ay(R) between massive E’s and Sa’s. The gradient in /.(* 
and (log age)L of SO’s is so similar to Sa’s galaxies, that they 
can result from the same formation process. 

- It is the Hubble type, not M*, that drives differences in the 
galaxy averaged age, and radial age gradients. These results 
indicate that the SFH and their radial variations are modu¬ 
lated primarily by galaxy morphology, and only secondarily 
by the galaxy mass. This suggest that galaxies are morpho¬ 
logically quenched, and that the shutdown of star formation 
occurs outwards and earlier in galaxies with a large spheroid 
than in galaxies of later Hubble type. 

From the comparison of the results with the theoretical pre¬ 
dictions from cosmological simulations, we conclude: 

- Major mergers are likely the main process building the cen¬ 
tral regions of ETGs. The metallicity gradient within 1 HLR 
is shallow compared with the theoretical expectation if mi¬ 
nor mergers are relevant in the growth of the central core of 
E’s and SO’s. In our results there is no evidence either of an 
inversion of (log age)r toward older ages beyond 1 -2 HLR, 
or of a steepening of the metallicity if these galaxies were 
growing in size through minor dry mergers. Massive galaxies 
probably accreted massive satellites that were able to retain 
their metal rich gas against winds, producing flatter metal¬ 
licity gradients (Hirschmann et al. 2015). Alternatively, the 
flattening of the metallicity radial profile can result from the 
quenching of star formation. When this happens, the metal 
cycle stops and only stars of that last star formation event 
remain. 

- Through the negative metallicity gradients, spirals show evi¬ 
dence of growing inside-out. These gradients are flatter than 
the predictions by the classical chemical evolution mod¬ 
els (Chiappini et al. 2001; Molla & Diaz 2005, e.g), but 
are similar to those measured above the Galactic disk. The 
largest gradient happens in intermediate types and interme¬ 
diate galaxy mass, as predicted by the Molla & Diaz (2005) 
models. However, Sbc galaxies have a (log Z*) jW gradient 
similar to the predictions by RaDES simulations (Few et al. 
2012; Pilkington et al. 2012a). This indicates that the feed¬ 
back recipes used in these simulations are able to recover 
realistic galaxies with small bulges. However, a larger set 
of cosmological simulations is required, covering from early 
type Sa to late Sd, and a large range of galaxy mass (from 
10 9 to 10 11 Mo), in order to confirm the general trend with 
the Hubble type found in this work. 
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Thanks to the uniqueness of the CALIFA data in terms of 
spatial coverage and resolution, large sample spanning all mor¬ 
phological types, and homogeneity and quality of the spectral 
analysis, we are able to characterize the radial structure of the 
stellar population properties of galaxies in the local universe. 
The results show that the Hubble sequence is a useful scheme 
to organize galaxies by their spatially resolved stellar density, 
ages, and metallicity. However, stellar extinction cannot discrim¬ 
inate so well between the different Hubble types. Stellar mass, 
although responsible for setting the average stellar population 
properties in galaxies, it is less responsible of the quenching pro¬ 
cesses. Morphology is, however, more strongly connected with 
the shut down of the star formation activity in the bulges and 
disks of galaxies. 
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Appendix A: Improvements resulting from the new 
reduction pipeline 

The quality of the starlight fits to CALIFA version 1.3c spec¬ 
tra was assessed in Cid Fernandes et al. (2014) by averaging 
R a = 0 A - M a residual spectra of 107 galaxies (~ 10 5 spectra). 
Inspection of these residuals revealed low amplitude (a few %) 
but systematic features related to unmasked weak emission lines, 
SSP deficiencies and data calibration imperfections. This exer¬ 
cise needs updating now that the reduction pipeline has changed 
to version 1.5. 

Fig. A.l summarizes the results of this re-evaluation. The 
plots show stacked R,\ = 0 A - M A residual spectra, in units of 
the median flux in the 5635 + 45 A window. The top panel shows 
results for the nuclear extractions, while the middle and bottom 
ones are built using spectra from zones within radial distances 
R — 0 - 1 and 1-2 Half Light Radius (HLR, computed in the 
same wavelength range), respectively. Residuals are colored ac¬ 
cording with the Hubble type of the galaxies. When all galaxies 
are stacking, the residuals are colored according with the spatial 
zone extracted. These subdivisions are presented in order to get 
a sense of how the residuals relate to position within a galaxy 
and its Hubble type, two central aspects of this paper. 

No matter which panel one looks at, the improvement with 
respect to version 1.3c is evident to the eye when compared 
to figure 13 of Cid Fernandes et al. (2014). The broad trough 
around H (3 present in the 1.3c spectra, for instance, is much shal¬ 
lower now. In fact, it is confined to late types (compare blue and 
red lines in the lower panel in Fig. A.l), indicating that its origin 
is not only related to calibration, but also to the SSP spectra of 
young stellar populations (as previously reported by Cid Fernan¬ 
des et al. 2005 for SDSS data). Residuals are also visibly smaller 
towards the blue, including the Call K line, which is now well 
fitted whereas in version 1.3c a small systematic residual sub¬ 
sisted. 12 The humps around 5800 A, on the other hand, are still 
present in version 1.5, particularly noticeable for outer regions, 
indicating that further refinement of the sky subtraction are war¬ 
ranted. 

In short, the spectral fits have improved substantially with 
the new reduction pipeline. We attribute this to the updated sen¬ 
sitivity curve used in version 1.5. A more extended discussion 
of these and other aspects of the data reduction are presented in 
Garcfa-Benito et al. (2015). 

Despite these changes, the stellar population properties de¬ 
rived from the spectral fits did not change much in compari¬ 
son to those obtained for version 1.3c data. The most noticeable 
changes were in mean ages, which become 0.1 dex older, and 
extinction, which is now 0.2 mag smaller on a global average. 


Appendix B: Base experiments and uncertainties 
due to SSP templates 

To derive the stellar population properties of these 300 CAL¬ 
IFA galaxies we have fitted -253000 spectra with the GMe and 
CBe using the cluster Grid-CSIC at the Instituto de Astroffsica 
de Andalucfa and the cluster Alphacrucis at IAG-USP Sao Paulo. 
Examples of the quality of the spectral fits as a function of the 
Hubble type and radial distance are presented in Fig.A.l. 


12 Because of this improvement, our starlight fits now start at 3700 A, 
whereas in previous articles in this series only the A > 3800 A range 
was considered. 



Fig. B.l. The global stellar MZR for 300 CALIFA galaxies is shown 
as dots, color coded by the morphological type. It is derived using 
(log(upper panels) and (log Z^)^ LR (lower panels), and the total 
stellar mass, M*, obtained with the GMe SSP models (left panels) and 
CBe (right panels). A mass-binned smooth mean relation is shown as a 
solid black or grey-black line. The MZRs obtained for SDSS galaxies 
by Gallazzi et al. (2005) and Panter et al. (2008) are plotted as brown 
and magenta lines, respectively, with dashed brown indicating the 16 
and 84 percentiles of Gallazzi et al. (2005). 

Appendix B.l: Mass-metallicity relation 

Here we want to find out how well our metallicity definitions 
follow a MZR that guarantees that galaxies like the MW or 
Andromeda (log M*(M Q ) ~ 11) have solar metallicity at the 
disk, while LMC-SMC-like galaxies (logM*(M Q ) ~ 9) have ~ 
1 /4Z 0 . We do this with mass-weighted and luminosity-weighted 
definitions of eqs. 2 and 3, and with the two sets of SSP models 
(GMe or CBe). 

Similarly, the correlation between the galaxy averaged stel¬ 
lar metallicities and the metallicity measured at 1 HLR, and the 
MZR, guarantee that the metallicity radial profiles scale with 
the galaxy stellar mass. However, the MZR in Gonzalez Del¬ 
gado et al. (2014a) was derived using the galaxy averaged stel¬ 
lar metallicity instead of the metallicity measured at 1 HLR, 
and using the mass weighted definition of the metallicity and 
only the results with the base GMe. For these reasons, we de¬ 
rive the MZR that results from using the mass-weighted and also 
the light-weighted definition of the metallicity, and the GMe and 
CBe bases. 

Fig. B.l shows the correlation of M* and (log Z*) , l / / A> (upper 
panels), and (log Z Jr ) 3/, ' A ’ (lower panels) for the GMe (left pan¬ 
els) and CBe (right panels) SSP models. The mass-metallicity 
relation found by Panter et al. (2008) and Gallazzi et al. (2005) 
are the magenta and brown lines, respectively 13 . Note that in the 
four cases, the metallicities are well in the range of the dispersion 
given by Gallazzi et al. (2005) (brown dashed line represent the 

13 The Gallazzi et al. (2005) and Panter et al. (2008) relations have been 
shifted by 0.25 dex to the right to account for the difference in IMF 
between their results based on models with Chabrier IMF and ours ob¬ 
tained with Salpeter IMF in the two left panels. 
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Fig. A.l. Upper panel: Residual spectra averaged for all the spectra (black line), nuclei (grey line) and spectra belonging to zones that are between 
0-1 HLR (pink line) and between 1-2 HLR (yellow line)). Middle panel: residual spectra averaged for zones inner to 1 HLR and for Hubble type. 
Bottom panel: As in the middle panel but for spectra of zones located between 1 and 2 HLR. 


16th and 84th percentiles of their distribution). To compare the 
general trend of these values, we have derived a smoothed mass- 
binned relation, represented by a black solid or grey-black line. 
As expected from the global MZR derived in Gonzalez Delgado 
et al. (2014a), base GMe and (log Z*)m predict stellar metal- 
licities for MW and LMC-SMC with the expected values. But 
(log Z*)/, gives a MZR that predict higher metallicities. The op¬ 
posite happens for the MZR using the base CBe, that gives mass 
weighted metallicities higher on average than the SDSS metallic¬ 
ities, but the MZR with (log Z*)/. goes close to the Gallazzi et al. 
(2005) relation, and also predicts stellar metallicities for MW- 
Andromeda-like and LMC-SMC galaxies with the expected val¬ 
ues. 


In summary, (log Z*)m with GMe and (log Z+)/. with CBe 
provide a mass-metallicity relation similar to the SDSS MZR, 
and predict metallicities between -0.7 and -0.4 dex for galaxies 
with mass between ~10 y and 10 10 M 0 , the expected values for 
LMC and SMC galaxies like, and solar for MW-like galaxies. 


Appendix B.2: Variations over bases GMe and CBe 

Although the two sets of models are built with the same stel¬ 
lar libraries (MILES and granada), base GMe stops at 1.5Z Q , 
while CBe goes up to 2.5Z 0 . Because MILES is built with stars 
in the solar neighborhood, it does not contain stars as metal rich 
as 2.5Z Q , so CBe results at over solar metallicity should be inter¬ 
preted with care. On the other hand, the central parts of galaxies 
can be as metal rich as 2-3Z Q , and the base GMe may be too 
low to fit spectra of these regions, leading to saturation effects. 
To avoid these problems, we have also fitted the spectra with an¬ 
other two set of SSPs, that are identical to GMe and CBe, but 
where for each metallicity bin, two SSPs of age 16 and 18 Gyr 
are added to our "standard" bases. These extra bases, that we 
name GMd and CBd, allow galaxies older than the age of the 
universe if their bulges are very metal rich. Furthermore, results 
at very low metallicity also must be taking with care of. MILES 
contain only few stars of metallicity below 1/100 Z Q . For this 
reason, Vazdekis et al. (2010) provide a safe age range for each 
metallicity bin, being the models with logZ*(Z Q ) < -1.7 only 
valid between 10 and 18 Gyr. This safety margin is provided to 
avoid the cases when, due to age-metallicity degeneracy, these 
old metal poor models fit young metal rich populations, that may 
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(log Z t ) L (Z e ) GMe 


Fig. B.2. Correlation between SP obtained with the GMe (x-axis) and CBe (y-axis bottom and middle panels; and GMe (x-axis) and GMd (y-axis). 
The average difference between the property in the y-axis and x-axis is labeled in each panel as A, and the dispersion as (<x). 


happen if the base does not include SSP younger than 100 Myr. 
Our fits do not suffer this problem because bases GMe and CBe 
both have spectra of ages as young of 1 Myr. 


Appendix B.3: Global results and uncertainties associated to 
SSP models 

To evaluate to which extent the spectral synthesis results de¬ 
pend on the choice of SSP models, we now compare the global 
properties derived with bases GMe and CBe. Using our pipeline 
pycAsso we have obtained the radial distribution of the stellar 
population properties for each galaxy with a spatial sampling of 
0.1 HLR. Here we compare the stellar population properties of 
the 0.1 HLR radially sampled points, instead of comparing the 
results obtained from the individual 253418 fitted spectra. Fig 
B.2 shows the results for a total of 6000 points corresponding to 
a maximum of 20 radial points (from nucleus to 2 HLR) for each 
of the 300 galaxies analyzed in this work. The figure compares 
the results for base GMe in the x-axis, with CBe in the y-axis, 
in the bottom and middle panels. The upper panels compare the 
results of GMe with GMd. Each panel quotes the mean A and its 


standard deviation, where A = property)CBe) - property(GMe) 
or A = property (GMd) - property(GMe). 

GMe-based /revalues are higher than CBe by 0.27 dex on 
average, reflecting the different IMF used. Apart from this off¬ 
set, the two stellar mass surface density agree to within 0.08 
dex. Mean extinction is also in good agreement with a disper¬ 
sion of 0.05 mag. Ages are higher in GMe than CBe by 0.14 
dex for (log age)i and 0.08 dex for (log age) with disper¬ 
sion 0.18 dex and 0.12 dex, respectively. In someway, this re¬ 
sult is expected since base GMe also differs from CBe in IMF 
and isochrones. The differences in opacities in the equation of 
state between Padova 2000 (GMe) and 1994 (CBe) tracks pro¬ 
duce somewhat warmer stars in the red giant branch in the for¬ 
mer. Thus, older ages are expected with GMe than with CBe. 
However, the metallicities are lower in GMe than in CBe by 
0.13 dex for (log (Z*/Z g ))m and very similar (on average) for 
(log (Z/Zq))l. In both cases, the dispersion is similar, 0.11 and 
0.13 dex, respectively. Note that for Z* > Z Q , the metallicities 
(weighted in light or in mass) are always higher with CBe than 
GMe, reflecting the saturation effects in the base GMe due its 
limitation to Z < 1.5Z 0 . The shift at under-solar metallicities 
may be reflecting the age-metallicity degeneracy, CBe giving 
higher metallicity and younger the ages. 
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The upper panels compare the results of GMe with the GMd. 
Here, we see two relevant effects. The results of GMd also differ 
from GMe in the the range of extinctions allowed to starlight. 
While with GMe and CBe starlight always assumes Ay > 0, 
with GMd , starlight can bluer the SSP spectra by up to Ay = 
-0.5 mag. This is allowed to avoid the effect of saturation at Ay 
= 0. The global effect is that ages can be 0.11 dex younger with 
GMe than GMd. Metallicity is not affected by this choice of Ay. 
But the extension to ages older than the age of the Universe in 
GMd has some effect on the metallicity above Z 0 , so metallicities 
are slightly lower and ages slightly older. 


zone; this is not possible for A,, (Fig. C.4), (log age)i (Fig. C.2), 
or (log Z t ) m (Fig. C.3), because these are intensive properties. 
Fig. C.4 shows how Ay changes across the CMD: the most lumi¬ 
nous galaxies are little affected by extinction, while Ay is higher 
in spirals of intermediate type and with blue colors. The mean 
(dispersion) Ay values at the nuclei, 1 HLR, and 2 HLR are 0.47 
(0.37), 0.19 (0.16), and 0.13 (0.13), respectively. 2D maps and 
the difference in the mean values at different distance indicate 
that stellar extinction shows radial gradients. 


Appendix C: 2D maps and tables 

Appendix C. 1: 2D maps of log 

Fig. C. 1 shows the M r vs. u-r CMD for the 300 CALIFA galax¬ 
ies of our sample. Each galaxy is represented by its 2D map of 
the log /.(* located at the position of its integrated M r and u-r 
values. In this plot the SFH is compressed into the present-day 
stellar mass surface density, that measures the end product of the 
SFH. Because our analysis accounts for extinction, these log /r* 
values and their radial variations are free from extinction ef¬ 
fects. Fig. C.l shows 14 clearly that log/r* correlates with M r , 
and spheroids are significantly denser than late type galaxies by 
one to two orders of magnitude at the center, and by one order 
of magnitude at distances 1-2 HLR. At the center, 2.0 < log/r* 
(M 0 pc' 2 ) < 4.7, while < log/i* (M 0 pc' 2 ) < 3.4 at 1 HLR, and 
1.0 < log a/* (M 0 pc' 2 ) < 2.9 at 2 HLR. 

Appendix C.2: 2D maps of (log age) L 

Similarly to Lig. C.l, Lig. C.2 shows the 2D maps of (log age)i. 
It portrays the correlation between the average age of the stellar 
populations and M r and colors, with the most luminous and red 
galaxies being older, while the bluest galaxies are the youngest. 
Gradients of the stellar population ages are also clearly detected 
within each galaxy in these 2D maps, and more remarkably 
in galaxies located in the green valley. At the center, 7.3 < 
(log age)L (yr) < 10.1, while 8.3 < (log age)/, (yr) < 10.1 at 1 
HLR, and 7.5 < (log age) L < 9.9 at 2 HLR. 

Appendix C.3: 2D maps of(logZ*) M 

Similarly to Lig. C.l, Lig. C.3 presents the 2D maps of 
(log Z*)m- It clearly shows the correlation between the stellar 
metallicity and galaxy luminosity, equivalent to the mass metal¬ 
licity relation. Gradients of the stellar metallicities are more 
clearly seen in these 2D maps in galaxies with intermediate lu¬ 
minosity (-22 < M r < -20). The most luminous galaxies have 
solar or over solar metallicity producing a visual saturation in the 
2D maps. The stellar metallicities range from (log Z*)m= -1-4 to 
0 . 22 . 

Appendix C.4: 2D maps of A v 

Similarly to Lig. C.l, Lig. C.4 presents 2D maps of Ay. Effects 
of spatial binning are visible in the Ay maps, where all the pixels 
within a Voronoi zone have the same value. These effects are not 
noticeable in the/i* images (Lig. C.l) because /r* is an extensive 
property, and the zoning effect was softened by scaling the value 
at each pixel by its fractional contribution to the total flux in the 

14 These 2D maps are the results from the GMe SSP base. 
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Fig. C.l. 2D maps of stellar mass surface density, /j*. Each galaxy is placed in its location in the u — r vs. M, diagram, where color and magnitude 
correspond to its global values. The 2D maps are shown with North up and East to the left. 



Fig. C.2. As Fig. C.l, but for images of the luminosity weighted mean age, (log age) L . 


Article number, page 29 of 42 












A&A proofs: manuscript no. Califa5_metallicity_vl2_8Jun 



Fig. C.3. As Fig. C.l, but for images of the mass weighted mean metallicity, (logZ*)^. 



Fig. C.4. As Fig. C.l, but for images of the stellar extinction, Ay [mag] 
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